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FOREWORD 


This report contains the results of the Study of Radar Pulse Compression 
for High Resolution Satellite Altimetry awarded Technology Service Corporation 
under Contract No. NAS6-2241 by the National Aeronautics and Space Administration 
Wallops Station, Wallops Island, Virginia. The study was conducted by Technology 
Service Corporation under the direction of Mr. Fred Nathanson as Program Manager 
with Dr. Richard P. Dooley as Assistant Program Manager, 

Successful implementation of this effort was due primarily to Mr. William 
Townsend, NASA/Wallops Program Manager, who provided considerable guidance and 
direction during the course of this program, 

A major contributor to this study was Dr, Lowell Brooks, Senior Scientist 
of Washington Operations, who performed the analysis of improved range tracking 
algorithms. The concept of a maximum likelihood processor, which has a signif- 
icant impact on the results of this study^ was originally suggested by Dr, Peter 
Swerling, President, Technology Service Corporation, Other researchers who 
contributed to this effort included Mr, James Bucknam who performed much of the 
system design calculations and analysis. Dr. Peter Tong who performed the study 
of binary phase code with digital processing. Dr, Glen Gray for the analysis of 
the linear FM generation requirements, Mr. Alexander Mac Mullen who developed 
the system implementation, and Dr, August Rihaczek who acted in an advisory and 
review capacity during the course of this program. 
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AB STRACT . 

A study is made of pulse compression, techniques applicable to a 
satellite altimeter having a topographic resolution of + 10 ctn, A systematic 
design procedure is used to determine the system parameters. The performance 
of an optimum, maximum likelihood, processor is analysed in a supporting 
study and provides the basis for modifying the standard split-gate tracker to 
achieve improved performance. Bandwidth considerations lead to the reconmenda- 
tibn of a full deramp STRETCH pulse compression technique followed by an analog 
filter bank to separate range returns. The implementation of the recommended 
technique is examined in detail. 
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PART I 


SYSTEM DESIGN 

1,0 INTRODUCTION AND SUMMARY 


The design of a high resolution satellite altimeter is described in 
this part of the report. The resulting design achieves all specified performance 
requirements. These performance requirements are given in Table L.l and the 
system design is summarized in Table 1,2, 

Section 2 describes a systematic design procedure for determining the 
system parameters. This procedure clearly identifies the tradeoffs among alternate 
designs and as such jprovides a basis for the selection of a design which can be 
achieved in the most efficient and economical manner. The dependencies between 
the various system parameters and performance requirements are examined in detail. 
It is shown that the results of the improved altitude tracking algorithms investi- 
gation, Part II, have a significant impact on the selection of the nominal 
system parameters. The form of the optimum (maximum likelihood estimate) processor 
led to modifications of the simple split-gate tracker which enable the performance 
requirement of 10 cm resolution at 10 meter wave height to be achieved with 
onboard processing. 

Section 3 examines the types of pulse compression, considered for the 
satellite altimetry experiment. Utilizing the set of required system parameters, 
the feasibility of each technique is examined in detail. Bandwidth considerations 
led to the selection of a full deramp STRETCH followed by an analog filter bank 
to separate range returns as the recommended technique. The state-of-the-art 
in the generation of linear FM signals, an essential part of the selected STRETCH 
technique, is examined in detail. 
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Recommendations concerning the implementation of the selected pulse 
compression technique are given in Section 4. While the reflective array compressor 
(RAC) is the preferred method for the generation of the 360 MHz 2,8 jj.sec 
linear FM signal, procurement would be required from MIT Lincoln Lab since there 
are presently no established vendors of the RAC line, and if obtained from in- 
dustry this approach would involve some development risk. Instead, a configura- 
tion using a lower bandwidth (60 MHz) delay line followed by frequency multi- 
plication (x6) is recommended for the ’’baseline design”. While TSC would prefer 
to see the implementation using RAC, perpendicular diffraction grating delay 
line (PPDL) ,and conventional surface waves in that order, all approaches are 
capable of meeting the easier specification of the ’’baseline design”. An 
analysis of the accuracy with which the ramp and deramp linear FM signals 
must be generated provides a linearity requirement of 0,2% and peak allowable 
frequency deviations of < 25 KHz (for one cycle of variation across the pulse). 

The implementation of the analog filter band for range processing is examined 
in detail. 

Table 1,1 SYSTEM PERFORMANCE REQUIREMENTS 


I. 

Geodetic Accuracy 

50 cm 

II. 

Topographic Resolution 

10 cm rms (7 cm allocated to 
system error) 

Ill, 

Wave Height Range: 

1-10 m crest-to-trough 


Accuracy: 

25 % 

IV. 

Correlation between pulses 

< 1/e 

V. 

Oceanographic phenomena of 
interest (maximum spatial 
frequency) 

.25 Hz 
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Table 1.2 


design summary 


I . Or bit Parameters 

a) Height 

b) ■ Inclination 

c) Eccentricity 

II. R adar Parameters : 

a) Antenna Beamwidth 

b) Pointing Accuracy 

c) Antenna Gain Peak 

Average 

d) Peak Power 

e) System Losses (other than 
processing losses in pulse 
compressor) 

£) Noise Figure 

g) Frequency 

h) Uncompressed Pulse Width 

i) Pulse Bandwidth 

j) Compressed Pulse Width 

k) Compression ratio 

l) PRI’max (Uncorrelated returns) 

m) PRF 

n) S/N (Single Pulse) 

o) Ocean Backscatter Coefficier_t 

f ST ^ 


556 km 

90° retrograde 
.0064 maximum 

3° (24 inch dish) 

ag = 1/2° 

34.9 dB 
34.25 dB 
2 KW 

5 dB 

5.5 dB 

13.9 GHz 

2.8 ixsec 
360 MHz 

3.0 nsec 
1000: 1 

1.8 KHz 

> 1.4 KHz 
10 dB 
+6 dB 
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Table 1.2 


DESIGN SDMMAilY (Continued) 


p) "'Receiver Weighting -26 dB Modified Taylor 

q) 'Vulse compression processing loss .55 dB 

r) 'Main lobe broadening due to tapering 23% 

Included in the design but considered optional. 


III. 


T racker Configuration 

Type 

Tracks 

Early Gate Width 
Late Gate Width 
Gate Separation 
Bandwidth 


IV. 


Pulse Compression 


Type 

Range Processing 
Filter Bank 
Number of Filters 
Frequency Range 
Filter Bandwidth 
Output Data Form 


Time required for full sampling 
A/D Sampling Frequency 


Modified Split-Gate 

Quarter power point of leading edge 

3.0 nsec 
48 nsec 

> 70 nsec 

1.0 Hz 

Full Derarap STRETCH 
Analog Filter Bank 
Discrete Passive 
30 

9.2 to 20.8 MHz 

385 KHz (3 nsec resolution) 

Two TTL parallel words 

A. Range Bin Number, 5-bits 

B, Range Bin Amplitude, 6-bits 

450 microseconds, max. 

< 1 MHz 
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Table L.2 DESIGN SUMMARY (Continued) 


V, Linear FM Generation 


Type 


Surface Wave 


Bandwidth 60 MHz 

Multiplier Chain X6 

Pulse Length 2,8 )j.sec 

Linearity of FM < 0.2% 

Peak Frequency Deviation < 25 KHz 


(one cycle of variation across pulse) 


One device used for both transmit and receive. 


VI, 


Wave Height and Return Shape Processing 


Type 

Averaging Time 
Number of Samples 
Sample Interval 


Averaged Samples of Power 
Return 

.1 sec 

30 

3,0 nsec 
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2.0 


SYSTEM PARAMETERS 


In this section, a nominal set of system parameters are determined 
for achieving the performance requirements of the satellite altimeter. The 
approach to this task is a systematic design procedure which clearly identifies 
the tradeoffs among alternate designs and as such , provides a basis for the 
selection of a design which can be achieved in the most efficient and economical 
manner. Obviously, the systematic design procedure was not employed until the 
final stage of the selection process. In fact, a major portion of the effort 
involved a detailed examination of the dependencies between the various system 
parameters and performance requirements. 

These studies produced two results which have a significant impact 
on the selection of the nominal system parameters. 

First , the effect of wave height on system resolution (RMS tracking 
error) has been determined. Previous expressions for RMS tracking error have 
assumed a smooth sea surface. It was found that, for a given resolution and 
signal“to-noise ratio, going from the smooth sea case (say 1 meter wave height) 
to a 10 meter wave height resulted in a factor of 100 increase in required PRF, 
This result is quite significant since the performance requirement of 10 cm 
resolution at 10 meter wave height requires considerable improvement in tracker 
performance compared to that originally envisioned for the smooth sea case. 

Second , a study of improved range tracking algorithms has shown that 
the performance of a split-gate tracker could be improved considerably by widen- 
ing the width of the late gate and having the early gate positioned well below 
the half power point of the return signal. These modifications to the simple 
split-gate tracker were indicated after a detailed examination of the form of 
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the optimum (maximum likelihood estimate) processor. Without this result, the 
performance requirement of 10 cm resolution at 10 meter wave height could not 
be achieved with onboard processing. The recommended tracker is a "modified 
split-gate" that tracks the ^ power point of the return signal. The width of 
the early gate is t and that of the late gate is 16 t, where t is the compressed 
pulse width. 

The selected set of system parameters were presented in Table 1.2 
of Section 1,0, The remaining portions of this section provide the rationale 
for the selection of these parameters. Section 2,1 describes the utilization 
of the systematic design procedure while the detailed calculations are presented 
in Section 2.2. 

2. 1 Systematic Design Procedure 

In Section 2.2, the dependencies between the various system para- 
meters and performance requirements are described in detail. These efforts 
are necessary prerequisites for obtaining a system design, but there remains 
a need for systematically organizing the design procedure to clearly reveal 
whether or not a particular design has been achieved in the most efficient and 
economical manner. 

In the following, it is shown that the design need not be based on 
trial and error methods, but can be accomplished with a systematic procedure. 

The guideline for such a design is to attain the given performance specifica- 
tions while minimizing equipment complexity. This is, of course, what every 
designer is attempting to do. The objective here is to provide a method by 
which this can be done systematically, to yield a design where the selection 
of every parameter value is justifiable. Moreover, in tlie process, an under- 
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standing is obtained of the cost of improving the performance should the specifica- 
tions be changed, and of the cost of achieving some critical performance parameter. 

As described above, the process of system design can be viewed as a 
multivariate ,constrained minimization problem. The problem is multivariate since, 
in general, several system parameters must be determined by the design procedure, e.g. 
pulse length, compression ratio, tracker bandwidth, etc* The constraints of the 
problem are provided by the system performance specifications (e.g. tracking 
accuracy, resolution, etc). Finally, the quantity to be minimized is the equipment 
complexity. That is, the best choice of system parameters is the set of parameters 
which, first, meet all the design goals, and second, can be implemented more simply 
than any other set of parameters which meet the specifications. 

There are at least two basic problems in rigorously solving the above 
minimization problem. The first is that of quantifying system complexity. This 
is an extremely difficult task, and it is felt that without a major effort any 
quantification formula would be, at best, highly controversial and, at worst, 
useless. Therefore, the judgement as to the system complexity implied by a set 
of parameters must be left to a competent engineer. The resulting design pro- 
cedure cannot then be mathematically rigorous (possibly to its advantage). 

The second basic problem, and the one which is addressed by the system- 
atic design procedure is that of determining all combinations of system parameters 
which will satisfy the performance requirements. It is felt that if these ’’feasible" 
system solutions are presented in an orderly manner, then the final judgement as 
to which is simplest can be made fairly easily. 

The design process can be summarized in 5 steps. 

1) Define precisely which parameters are to be determined. 
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2) Define those input and system parameters which 
have previously been determined from other con- 
siderations. 

3) Determine the various dependencies between the 
system parameters, the parameters, and the 
system specifications. Make a precedence table 
to indicate these dependencies. 

4) Use a procedure developed by Steward [l], to re- 
order the precedence matrix and to develop a flow 
chart which shows the order in which the parameters 
are to be determined. This yields a systematic 
procedure for exhaustively examining all feasible 
sets of system parameters, 

5) From the flow chart developed in 4) , compute and 
present tables of feasible solutions, and then 
select the set of parameters which gives minimum 
complexity. 

2.1.1 Application to altimeter design 

Step 1 . As defined in the work statement [2], and from basic consider- 
ations, the parameters which must be determined to define the altimeter are given 


in Tab le 2.1,1. 


TABLE 2.1.1. Altimeter System Parameters to be Determined 


Min 


J 

V 

y 

■j 

y 


Max 

y 


y 


y 


Parameter 

Spatial frequency 
Number of pulses integrated 
Signal- to-noise ratio 
Compressed pulse length 

Compression ratio 
Tracker bandwidth 
Pulse rep. freq. 


Symbol 

SF 

N 

S/N 

T 

CR 

^T 

PRF 


In general, all of the parameters would have a range of values which 
lead to feasible system solutions. However, in many cases only one end of the 
range will have any impact on the system design. For example, consider the com- 
pressed pulse length. Although in principle, there may be a minimum pulse length 
which will meet the system specifications, this will have little impact 
on the design. That is, the complexity of the pulse compression system increases 
rapidly with system bandwidth. Thus the designer will always want to minimize 
the system bandwidth, or equivalently, use the longest compressed pulse he can» 
Thus he is only interested in the maximum pulse length which will still meet 
the system specifications. 

By similar arguments, it can be shown that- only the maximum spatial 
frequency and the minimum- number of pulses integrated, signal-to-noise ratio, 
compression ratio, and tracker bandwidth are of concern from a design standpoint. 
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In the case of PRF, the minimum PRP is of concern from a system com- 
plexity standpoint, however, since the return for a very high PRP becomes 
correlated, there is an upper limit on the PRF which must be considered. 

S tep 2 . Table 2.1,2 outlines the performance specif ications , and 
system parameters which have been previously determined from other considerations. 

S tep 3 , Table 2,1.3 shows the precedence matrix of interrelations 
between parameters. This matrix was obtained from the various dependencies 
outlined in Section 2.2. The dependencies between the parameters can be seen 
by reading down columns of the matrix, and "X” indicates a dependency. For 
example, the maximum spatial frequency to be tracked can be computed from tables 
of the oceanographic phenomena of interest (OPI) and the satellite orbit parameters 
(OP). 


Similarly the minimum compression ratio required (CR) can be found 
once the input radar parameters (R) and the sea surface crosssection (t*^) are 
given, and after the minimum signal-to-noise ratio has been determined. 

Step 4 . The reordered PTBD precedence matrix using Steward's algorithm 
is given as Table 2.1.4, The purpose of Steward's algorithm is to put the matrix 
into block upper triangular form. By reordering the matrix so that the "X's" 
fall above the diagonal, it becomes immediately apparent which parameters must be 
determined first. One can then develop a flow chart as in Fig, 2,1,1, which allows 
a systematic development of all feasible sets of parameters which satisfy the 
performance constraints. 

For example, the max spatial frequency depends only on the input para- 
meters OPI and OP (Table 2,1,3) and not on any other system parameters. Thus it is 
determined first. From SF and the wave height, the maximum pulse length and 
the bandwidth of the tracker are determined. Third, from T the PRF is determined. 

ul2.X 
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.‘ABLE 2.1,2 Input Parameters, and Systems Parameters 
Which Have Been Previously Determined. 


B. 


Requirements from Specification 

Value 

I 

Geodetic Accuracy 

GA 

50 cm 

II 

Topographic Resolution 

TR 

(7 cm 

10 cm rms 

allocated to system error) 

III 

Wave 

Height Range: 

WH 

1-10 m crest-to- trough 



Accuracy : 


25% 

IV 

Correlation between pulses 

P 

^ l/e 

V 

Oceanographic phenomena of interest 

OPI 

Table 2.2.1 

System Parameters which have been specified 



I 

Radar parameters; 

R 



a) 

Antenna Beamwidth 


3“ 


b) 

Pointing Accuracy 


” 1/2° 


c) 

Antenna Gain Peak 


34.9dB 



Average 


34.25dB 


d) 

Peak Power 


2 KW 


e) 

System losses (other than processing 
losses in pulse compressor) 


5 dB 


f) 

Noise Figure 


5.5 dB 


g) 

Frequency 


13.9 GHz 


h) 

*Pulse compression processing loss 


.55 dB 


i) 

*Main lobe broadening due to tapering 

23% 

II 

Ocean Backscatter Coefficient 

o 

n 

+6dB 

III 

Orbit Parameters 

OP 



a) 

Height 


556 km 


b) 

Inc lination 


90^ retrograde 


c) 

Eccentricity 


,0064 maximum 

* F 

or an assumed ’iodified Taylor weighting, -2! 

>.7 aB peak 

side lobe • 
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TABLE 2.1.3 Precedence Matrix 






Parameters 

To Be 

Determined . (PTBD) 





SF ^ 
max 

N . 

min 

S/N . 
min 

Xiax 

CR B 

min T niin 

PRF . 
min 

PRF 

max 

p 

GA 








A 

TR 


X 

X 





R 

A 

WH 


X 

X 

X 




M 

i 

P 







X 

E 

DPI 

X 



Weak 




T 









E 

R 

R 





X 



S 

o 

a 





X 




OP 

X 



Weak 



X 


SF 

//// 



Weak 

X 

X 



1 N 


//// 

X 



X 



1 S/N 


X 

//// 


X 




i 

T 

1 




//// 

X 


X 


1 

1 CR 





//// 









//// 

X 



PRF 

tnin 






//// 



PRF 

max 







//// 


TABLE 2.1.4 

Re-Order Precedence Matrix 





max 

T . PRP 

max T min max 

N . 
min 

S/N . 
min 

CR . 
min 

PRF . 
min 


Now, in Table 2. 1.4, the number of pulses integrated and signal-to-noise 
ratio occur as a block on the diagonal. Thus, these parameters must be varied 
simultaneously since they cannot be factored into a precedence order. Having 
chosen them, the compression ratio is determined next, and the minimum 
FRF is determined last. 

Step 5 . From the performance criteria, the first four parameters 

(SF, T, B , FRF ) are determined almost uniquely. They are given in Table 2.1.5. 

X TU^X 

As noted in that table, the max spatial frequency is determined by the Gulf Stream, 
and is ,75 Hz for an extreme case. The pulse length is determined primarily 
by the minimum wave height resolution criterion. 

Thus in order to measure wave height down to 1 m, the pulse length 
must be no greater than about .5 m (3 nsec). In order to keep tracking biases 
down, the min bandwidth of the tracker is made 4 times the max spatial frequency. 
The max FRF is determined from t to be 1,8 KHz, 

The remaining parameters are presented in tabular form since they 
must be varied simultaneously. Three tables are presented which correspond to 
three types of epoch tracking systems, i.e. thestandard split-gate tracker, the 
modified split-gate tracker, and the MLE tracker. These systems are ranked in 
order of increasing complexity. 

The system selected based on considerations of system complexity and 
development risk is boxed in on Table 2.1,7. 
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TABLE 2.1.5 Parameters Shown to Have Nearly Unique Values 


Parameter 


Selected Values 


max SF 


max T 
min 


.75 Hz Extreme Gulf Stream 
•1 “.25 Hz Typical Gulf Stream 

3 nsec 

3 Hz Extreme Gulf Stream 
1 Hz Typical Gulf Stream 


max PRF 


1.8 KHz 


TABLE 2.1,6 Half -Power Split-Gate Tracker 



T ^ = 16t 
gL 




WH = 10 meters 

k = .5 (tracks half -power point) 


0=7 cm. 
T 

Continuous model 










PRF . 
mxn 

S/N (dB) 


N 


CR 

<3 6^=1 

0 

19-4 

* 

10^ 

96 

19.4 KHz 

5 

6.0 

* 

103 

303 

6.0 

10 

3.4 

* 

10^ 

957 

3.4 

15 

2.7 

* 

10^ 

3030 

2.7 

20 

2,5 


10^ 

9570 

2.5 


16 



Table 2.1.7 Modified Split-Gate Tracker 


T = T , T = 16T WH = 10 meters 

gE gL 


k = 

.25 (tracks quarter-power 

point) 

= 7 cm. 

T 


Continuous 

model 


PRF . 
min 


S/N 

(dB) N 

CR 

(B = 1 Hz 


0 

13.5 * 10^ 

96 

13.5 KHz 


5 

3.0 * 10^ 

303 

3,0 


r 

1 10 
L 

1.4 * 10^ 

957 

1.4 

1 

I 

15 

1.0 * 10^ 

3030 

1.0 


20 

0.9 * 10^ 

9570 

0.9 


Table 2,1.8 

MLE Tracker 0 " 

max 

16, WH - 

10 meters, cr^ = 7 

T 

cm. 




PRF . 
min 

PRF . 
min 

S/N 

N 

CR 

(a = 1 Hz 

T 

(a B^ = 3 Hz 

0 

1.9 * 10^ 

96 

1900 

5 700 

! 3 

600 

303 

600 

1800 1 

‘ 10 
1 

345 

957 

345 

1 

1050 ' 

J 

15 

270 

3030 

270 

810 

20 

240 

9570 

240 

720 
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A system based on the half-power split-gate tracker is preferred since it 
is easiest to implement and its characteristics are well understood from the Geos-C 
program. However, as an examination of Table 2.1.6 shows, a system based on 
the half-power split-gate tracker does not meet the performance specification 
unless it operates at PRP’s greater than 2,5 kHz. But at this rate, the maximum 
PRF of 1.8 kHz is exceeded, and the pulses become correlated. Therefore, it is 
unlikely that the half-power split-gate tracker will meet the specifications. 

The modified split-gate tracker brings the PRF down to an acceptable 
level for S/N greater than 10 dB while the MLE has acceptable PRF's at all S/N 
greater than about 5 dB, 

Note that in both cases, as the S/N increases, the required compression 
ratio increases rapidly. For compression ratios greater than about 1000, the 
pulse compression unit becomes a higher risk development item. Thus the sets of 
feasible solutions are reduced to the portions of Tables 2.1.7 and 2.1.8 
corresponding to a modified split-gate tracker operating at about S/N = 10 dB, 
and a MLE operating in the range of S/N =5 to 10 dB. Of the two, the MLE places 
less stringent requirements of the transmitter duty cycle (due to the low PRF), 
However, from a development standpoint, the higher performance transmitter is 
thought to be a lower risk item since one which meets the requirements [3] is known to 
exist. The MLE, however, must be considered as high risk since only theoretical 
performance calculations have been made, and no development work has been done. 

The system chosen is summarized in Table 2.1.9, 

The epoch tracker is a modified split-gate tracker which tracks the 
quarter power point of the return, the early gate width is 3 nsec, and the late 
gate width is 48 nsec. 
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TABLE 2,1,9 System Parameters-Determined from The Design Procedure 


Parameter 


Symbol 


Value 


Max Spatial Freq. 

Min Number of Pulses 

Min Signal-to-Noise 
Ratio 

Max, Compressed Pulse 
Length 

Min Tracker Bandwidth 
Min Pulse Rep, Freq, 


Max Pulse Rep, freq. 


SF 

N 

S/N 


T 


PRF . 
mxn 


PRF 


max 


.25 Hz 
1400 

10 dB 

3 nsec 

1 Hz 
1400 Hz 

1800 Hz 
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2,2 


Parameter Calculations 


The dependencies between the various system parameters and performance 
requirements presented in the previous section are now examined in detail. In 
addition to the design equation or rationale used to calculate the parameter values 
presented in Table 2,1.9, details concerning the specification of the antenna 
parameters and receiver weighting are also included for completeness, 

2.2.1 Survey of oceanographic and geodetic signals of interest 

A survey was made of those oceanographically and geodetically induced 
variations in satellite altitude which the altimeter should be designed to track. 
The survey determined for each such variation, the characteristic amplitude, 
rise time, and maximum rate of change of altitude. The results of the survey are 
shown in Table 2,2.1, The briefest rise time (1.3 - 6,5 sec) would be caused by 
such phenomena as boundary currents and eddies (e.g,, the Gulf Stream) and higher 
frequency undulations of the geoid, while the maximum rate of change of altitude 
that could be expected would be due to the eccentricity of the orbit itself 
(f^ + 50 m/ sec) . 

2.2.2 Tracker bandwidth 

The tracker bandwidth should be sufficiently wide such that several 
uncorrelated tracker outputs are obtained during the shortest rise time in Table 
2,2,1, A bandwidth of 1 Hz would meet the criterion at all but the worst case 
Gulf Stream (10 km width stream, perpendicular intersection of stream and orbit). 

A bandwidth of 3 Hz, while satisfying the worst case Gulf Stream, would result in 
an excessive PRP (4.2 kHz) for the \ power split-gate tracker. Consequently a 
tracker bandwidth of 

\ 1 

is chosen. 
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TABLE 2,2.1 Survey of Geodetic and Oceanographic Signals of Interest 


Phenomenon 

Spatial 
Exfpnt (km) 

Amplitude 

(meters) 

Max Range 

iV 

Rate (m/sec) 

Rise 

* 

Time (sec) 

Western Boundary Currents 





Typical Gulf Stream 

100 

1.0 

.08 

13 

Worst Case Gulf Stream 

10-50 

1.0 

.15-. 8 

1.3-6. 5 

Boundary Current Eddies 

100 (near stream) 

.35 

.03 

13 


200 (far away) 

.05^ 

.002 

25 

Open Ocean Currents 

500-1000 

.10 

.0008-. 0015 

65-130 

Coastal Sea Level Slope 

2200 

.60 ■ 

.002 

300 

Difference in Sea Level 
(East/West) 

— 

.60 

— 

■ 

Tsunamis 

50 (open seas) 

• 

U) 

. o 

1 

. 

o 

.05-. 08 

6.5 

Geoid Undulations (such 
as Puerto Rican and 
Venezuelan Trenches)' 

100-150 

10-20 

.5-1.5 

13-20 

Waves (sea and swell) 

7,6 km grid 
spacing at 

1-10 (peak 
to trough) 



’ 


1 sample/sec 




Orbit Eccentricity 

one r evo lu t i on 

+ e(r + h) 
— e 

= + 44.4 

48.9 

2800 

(orbital 

half- 

period) 


* 

at h 


300 n.mi,, = 7,63 km/ sec 
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2.2.3 C ompressed pulse width 

The compressed pulse width ( after any tapering effects) is chosen such 
that the minimum significant wave height to be measured (1.0 meter peak- to- trough) 
is at least twice the compressed pulse length; i.e,.^ 2 samples in the rise time 
of the sea echo leading edge. While 3 samples might seem more desirable, the 
resulting bandwidth (~500 MHz) is considered excessive. Thus, a compressed pulse 
length of 

T = 3 nsec (after tapering) 
is selected, 

2.2.4 Pulse return decorrelation time 

In order to determine the maximum PRF such that pulse- to-pulse 
fluctuations in the sea echo are uncorrelated, the decorrelation time of these 
fluctuations must be found. This decorrelation time will be determined by three 
effects: 

1) The Doppler spreading of the spectrum of the compressed return 
pulse envelope due to the horizontal velocity of the satellite; 

2) The Doppler spreading due to the random velocities of the 
scatterers (wave spray) ; 

3) The degree of overlap of the footprints associated with 
successive pulses. 

For the orbital parameters of interest, and for a compressed pulse width of 3 
nanoseconds, the first effect, Doppler spreading due to horizontal satellite 
velocity, is the dominant effect. 
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The doppler spreading is a function of the effective footprint size, 
which in turn is a result of antenna shaping, surface shaping, and pulse shaping 
functions. For a satellite altitude of 300 nautical miles, antenna beamwidths 
of a few degrees, and pulse lengths of a few. nanoseconds , the footprint is 
pulse- limited , as shown in Figure 2,2,1, and hence the decorrelation time is a 
function of the compressed pulse length. 


Assuming uniform return from the pulse-limited footprint, the doppler 
spectrum is well approximated by a uniform power spectrum between + f^ as shown 
in Figure 2,2.2. Then the correlation function of the square-law envelope detected 
output is given by 


R;^(T) = 


sin^ (2n f T) 

Q 

(2tt f T)^ 
o 


with the first zero occurring at 


T =J— = ^ 

1 2 f ^ V, sin 0i 


where 




half the pulse- limited beamwidth 


= cos 


-1 


h + ^ 


Assuming uniform vertical wave-motion Doppler from + 3 m/sec, the 
square law detector output correlation due to vertical wave-motion is simply 
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Fig. 2.2.2D0F1LER SFECTRUM 'OF RETURN PULSE ENVELOPE DUE TO 

HORIZONTAL VELOCITY OF SATELLITE 


h = 300 nautical miles 
V =7,63 km/sec 




The correlation proportional to percentage overlap of the footprint is 



For a 3 nanosecond pulse, mean satellite altitude of h = 300 nautical 
miles, and = 7,63 km/ sec, these decorrelation times are 

= ,55 msec 
= 1,8 msec 


- 185 msec 
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2.2.5 


Maximum PRF 


The maximum PRF to ensure uncorrelated pulse- to-puise fluctuations is 
the inverse of the decorrelation time determined in the previous section. Ignor 
ing effects of wave spray and pulse overlap, the maximum PRF is given by 


PRF = — 
max T. 


=1.8 kHz 


2.2.6 S/N and N 

The relationship between the number of pulses required, N, and the 

single-pulse signal-to-noise ratio, S/N, at the output .of the. pulse compressor is 

determined by the maximum allowable random error in the tracker o^, the maximum 

significant wave height WH (peak- to- trough) at which this accuracy must be 

nid.x 

achieved, and the configuration of the tracker. The general form of this rela- 
tionship is 


o ViT 

f(S/N) 

where f(«) is some non-linear function determined by the tracker configuration, 
and T is the rise time (expressed in meters) of a linear fit to the leading edge 
of the sea echo. T is given by 


(3,1) WH 


T = 


max 
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Figure 2.2.3 shows a plot of f(3/N) for a quarter-power split-gate tracker, with 

an early gate width matched to the compressed pulse length and a late-to-early- 

gate width ratio of 16. For a S/N of 10 dB, f(10 dB) = .34. For WH = 10 meters, 

max ’ 


and cm, the number of pulses, N, is 


N = 


.2 


f (S/N) 


r 


0 . 1 ) ( 10 ) (. 34 ) 

(4) (.D//2 


1400 


2.2.7 Required compression ratio 

The compression ratio ,CR, is chosen to provide the desired pulse com- 
pressor output signal- to-noise ratio. It can be computed by the radar range 
equation, in which the sea surface backscatter is accounted for by a point target 
with cross-section equal to the pulse limited footprint area times the backscatter" 
coef f icient ,a° . 

For this model, the pulse-compressor output signal-to-noise ratio 

is given by 


where 


(S/N) IF 


P,(G^ V' ^ 

3 4 

(4n) h kT F 

o n 


= peak transmiter power 
= boresight antenna gain 
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in (meters) for MLE and Split Gate T-^a^ke^s (C^ “2.5 meter, d ■ .5 meter) 




a 


L 

s 

h 

F 


IF 


gain loss associated with pointing errors 

wavelength 

target cross-section 

system losses 

additional losses due to tapering on receive only 

altitude 
noise figure 

IF bandwidth 


The target cross-section is given by 


Cf 


O' 

a 


TT chT 


where 

CJ^ backscatter coefficient 
c = speed of light 

T = compressed pulse length (after tapering) 

and ncliT is the area of the pulse-limited footprint. 

The IF bandwidth, in terms of the after-taper compressed pulse length 

T , is 


B 


IF 


1 + O' 
T 


where ot represents the main- lobe broadening over uniform weighting. 
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Solving for CR 


CR = 


4(4n)^ kT (l + cx) F 

o n 

Ft (G^ a° C t2 


(S/N) 


= C (S/N)^^^ 

C is evaluated as follows; 

4(4tt)^ = 28.0 dB 

= (556 X 10^)^ = 172.4 dB 
-21 

kT = 4 X 10 watt- sec = -204 dB watt-sec 

o 

F = + 5.5 dB 

n 

1 + Of = 1 • 23 — .9 dB 
= 2kw = 33 dBw 

=2 (34.9) dB = 69.8 dB 

=2 (-.65) dB = -1,3 

= (.02157 m)^ = -33.3 dB 
= + 6 dB 

C = 84.8 dB m/sec 

L = - 5 dB 

s 

= -.55 dB 

■ ; t 

= (3 nsec)^ = -170.5 dB sec^ 

Thus, 

C = 19.9 dB, 

and the required CR for a 10 dB is 960. 

- 31 - 


out 



2.2.8 


R equired PRF 


The PRF required to average N pulses is determined by the tracker 

bandwidth, 


PRF = N • 

req'd L 


For = 1 Hz, N = 1400, 


PRF , , = 1,4 kHz 
req'd 


Note: 


= 3Hz corresponds to = 4.2 kHz which exceeds the maximum 

PRF for uncorrelated returns. 


2.2.9 Receiver weighting 

The effect of range sidelobes on altimetry bias and wave height 
measurement has been examined in [4j, These results, extended to more general 
cases and corrected for a computational error, are summarized in Fig, 2.2.4. 

There, both waveform and tracker bias versus RMS wave height data (normalized to 
the compressed pulse width, t^) are presented for uniform and 25 dB modified 
Taylor receiver weighting. The tracker used for these computations was a 
standard % power split-gate with early and late gate widths both matched to the 
compressed pulse width. For direct comparison, the data for the mean power 
response biases with and without receiver weighting are also presented in Tables 
2.2,2 and 2,2,3, respectively. As shown, there is no appreciable change in bias 
with or without receiver weighting and the bias that does arise from these range 
sidelobes is quite small. While computations have not been made for the "modified 
Split-Gate", quarter power tracker recommended in Section 2.1, it is felt that 
these biases (while somewhat larger) would remain less than 1cm at = 2.5 meters. 
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Figure 2.2*4 


Normalized bias vs waveheight due to 
range sldclobcs only . 


TRACKER BIAS-UNIFORM 






TABLE 2,2.2 Mean Power Response Bias -25 dB Modified Taylor Weighting 



10 

5 

4 

2 


1 

.25 

.065 

.036 

.033'‘ 

.025* 


.5 

.073 

.061 

.051'' 

a. 

.046'' 

[biasj in 

1.0 

.105 

.093 

.093* 

.089 

centimeters 

2.5 

.233'" 

.225 

.228'' 

.227 



'ABLE 2,2.3 Mean Power Response Bias - Uniform Weighting 


N. T 

\^(nsec) 

10 

5 

4 

2 


.25 

1 

.079 

.049 

. * 
.043 

* 

.035 


.5 

.093 

.069 

.069'' 

.060* 

|bias| in 

1.0 

.142 

.122 

* 

.120 

.110 

centimeters 

2.5 


.292 

. 300'' 

.292 



* 


Graphical Interpolation 



A word of caution- the biases described. above are only those due to 


range sidelobes causing the mean power return to differ from the ideal impulse 
response; i.e,, differ from an asymmetrical function at the half power point. 

In fact, as shown in Part II, Section 2, tracker bias is, in general, a function 
of wave height and signal- to-noise ratio for an ideal impulse response. 


If ,as in GEOS, the average voltage on the No, 8 waveform sampler were 
used as a measure of wave height, the response , sensitivity to wave height shown 
in Fig. 2,2.5 results. Here, for a compressed pulse width of 10 nsec, the slope 
of the weighted and unweighted response are essentially the same. For a 3 nsec 
compressed pulse width the difference would be even less. Thus the only degrade 
tion in performance, caused by receiver weighting, would be due to the usual re- 
duction in S/N ratio (about 207o for the 25/ dB Modified Taylor). 

In summary then, it would seert that the 13 dB sidelobes associated 
with no receiver weighting cause no problems as far as bias or wave height 
measurement are concerned. On the other hand , a limited amount of receiver 
weighting (say 20 or 25 dB sidelobes) provides a slight reduction in bias at 
the expense of a slight decrease in S/N. Smaller sidelobes have not been con- 
sidered since phase errors and other tolerance problems associated with the 
physical realization of a pulse compression technique can (and often do) cause 
the far out sidelobes to be much larger than the designed value when more than 
25 dB reduction is attempted. This being the case, there just doesn’t seem to 
be any good reason for either recommending or rejecting receiver weighting. As 
such, a 25 dB Modified Taylor receiver weighting has been included in the design 
but should be considered optional. 
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Figure 2.2,5 

AVERAGE SAMPLE VOLTAGE ON NO. 8 
WAVEFORM SAMPLER VERSUS WAVEHEIGHT 






2.2.10 


Antenna parameters 


The antenna gain specified in Table 2,1.2 cannot be obtained with an 
18” dish at 65% efficiency. Although this is about the highest efficiency which 
can practically be achieved with a parabolic dish, a higher gain at a given beam- 
width can be achieved by using a larger dish (smaller f/D) with effectively a 
heavy illumination taper to give the desired beamwidth. The efficiency of the 
larger dish will be even lower (e,g.«55%), but the gain will approach that of a 
uniformly illuminated dish of the same beamwidth. Using this approach, it is 
possible to realize an 87%. "efficiency" relative to the area of a uniformly 
illuminated dish of equal beamwidth, with off-the-shelf antenna. Thus, a realiz- 
able antenna gain at a 3° beamwidth would be 34.9 dB. 

However, the gain must also be corrected for losses due to beam point- 

2 Y 

ing errors. These losses can be accounted for by replacing with G , the average 
value of two-way gain. This average value can be computed by assuming a gauss ian 
beamshape and gaussian distributed pointing errors. Specifically, 

= J G^^G^ce, 0) P (e, 0) d0d0 

where 

= boresight gain (one-way) 

G(0, 0) = normalized antenna pattern (one-way) 

^-rr(0^ + 0^) /b2 

B = equivalent beamwidth 
p(Q> 0) - probability density of pointing errors 
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1 -(02 + 02 )/ 2 a 2 

= 2 ^ 

2 to 

~2 

Substituting into the expression for G yields 

2 9 a2/ 2 2 

= g/ {;e-2"0 ,3 } 

/2m 

2 2 1 
= (1 + 4na )“J- 

' 

= v" 

Thus the loss due to pointing errors is 
L„ = (1 + 4"a )*^ 

For 2a = 1^, this loss is -.65 dB at B = 3*^, 
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3 . 0 SELECTION OF PULSE COMPRESSION TECHNIQUE 

In this section the types of pulse compression that were considered 
for the satellite altimetry experiment are examined in detail. Utilizing the set 
of nominal system parameters from the previous section, each type of pulse com- 
pression is considered on the basis of feasibility, complexity, efficiency and 
stability. Bandwidth considerations led to the selection of - a full deramp 
STRETCH (similar to ALCOR) followed by an analog filter bank to separate range 
returns - as the recommended technique. 

3.1 Summary of Candidates 

The following pulse compression techniques have been considered 
for the satellite altimeter: 

1) Binary Phase Coding 

2) Linear Frequency Modulation 

3) Hybrid (analog/digital) 

4) STRETCH-ALCOR 

The binary (or polyphase) coding techniques are described in Section 
3,2, There it is shown that while the digital techniques have the desirable property 
of being able to change waveform (compression ratio), wave height data cannot be 
obtained with a simple 1 bit I, 1 bit Q system but requires a "multi-bit” decoder. 

For the required bandwidth, the complexity of even a 2 bit I, 2 bit Q (plus sign) 
system is considered to be pushing the state-of-the-art beyond 1973 technology and 
thus, the phase coded technique is not recommended at this time. 
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Section 3.3 considers the linear FM technique which is certainly the 
simplest and most widely used form of pulse compression. The problem with this 
method is found to be not so much pulse compression per se but the digital read- 
out of the resulting 300-360 MHz signal. . Even with sample-and-hold circuits, 
digitizing a 300 MHz signal with 6-8 bits per word is not considered practical 
and the linear FM (full compression) technique is not recommended. Section 3,3 
also contains considerable material on the state-of-the-art for the various 
methods of generating a linear FM signal since these signals are an essential part 
of the more general STRETCH. ALCOR configuration. 

Several hybrid analog/digital techniques are examined in Section 3.4. 
The hybrid of Barker code and linear FM is used to illustrate the fact that such 
techniques, while useful for increasing achievable compression ratio, in general 
require processing at the full signal bandwidth and hence are not recommended. The 
use of a binary phase coded waveform in a tracking mode, "cross-correlator", is 
shown to be capable of performing the altimetry but provides little or no informa- 
tion from which wave height can be accurately determined. 

The STRETCH-ALCOR techniques are examined in Section 3.5, The general 
technique is shown to be capable of reducing the bandwidth of the compressed signal 
and hence the A/D conversion requirement. Bandwidth and delay requirements for the 
STRETCH dispersive line and sampling frequency for the A/D convertor are given as 
a function of STRETCH ratio (SR). The fall deramp (SR = oo) followed by an analog 
filter bank to separate range returns is recommended over partial deramp (SR > 1) 
since this technique requires the least sampling frequency (< 1 MHz) for the A/D 
convertor and a single dispersive line for the generation of both the transmit and 
receive linear FM waveform. Digital filtering is not recommended since the A/D 
conversion would require a 21,4 MHz sampling frequency as compared with < 1 MHz 
for the analog filter bank. 
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3.2 


Binary Phase Coding 


The applicability of binary (or poljrphase) coding to satellite alti- 
meters depends upon three factors; 

1. The elimination of range-doppler ambiguities inherent 
in a linear FM or Chirp waveform, 

2. The availability of digital microelectronic signal 
processing techniques to directly give digital 
information on altitude and "sea state", 

3. Flexibility to change waveform with digital 
implementation. 

These factors can be discussed separately. The binary waveform is 
often chosen when the velocity of the vehicle or target is so uncertain that an 
absolute determination of time delay is impossible without an absolute deter- 
mination of relative radial velocity. The error in time delay (At^) measure- 
ment is proportional to the ratio of the doppler uncertainty (Af.) to the FM 

d 

dispersion of the waveform (AF) , 




AF 


T 


where T is the time dispersion. 

The source of doppler error is the result of the uncertainty in 
the eccentricity of the orbit. Since the maximum f^ is specified at 5 KHz, 

O 

the uncertainty should be of the order of 0,5 x 10 Hz, If we assume AF 300 
MHz and T ks 3 x 10 sec, then At^ % ,6 x 10 sec. Thus the ambiguities in 
the linear FM waveform do not seem to cause a problem and the choice of wave- 
form depends on ease of implementation. 
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The .binary waveform can be implemented, either in analog or digital 
form. With thi^ pJarameters of interest, an analog implementation would probably 
also use surface wave techniques. This is illustrated in Figs. 3.2-1 and 3.2-2. 
However, for a given time- bandwidth product, it is somewhat more difficult to 
-implement the .binary phase waveform with surface wave devices [2], [3] and [4], 
Since the potential advantage of the binary technique is that the output could 
be directly in digital form, there seems to be little value in further discussion 
of analog techniques. 

The simplest and most convenient form of decoding for a binary phase- 
cqded waveform is the "one-bit I, one-bit Q" system shown on Fig. 3.2-3. (From 
[l]) the received signal is mixed with the transmit carrier and only the polarity 

f ' ■«. 

of the bipolar in-phase and quadrature signals are entered into high speed shift 
registers, digital comparators and adders, (These must all work at a clock rate 
equal to the bandwidth of the transmitter waveform). Using digital adders, the 
maximum output is equal to the time -bandwidth product in each channel for a 
single point target. This assumes that the number of stages is equal to the 
time bandwidth product. Many decoders of this nature have been built with 5-20 
MHz bandwidths, and a few experimental models with higher bandwidths have been 
constructed. It seems possible to get to over 120 MHz bandwidth with MECL 

circuits, but there is question as to the practicality at 300 MHz. This is 

explored further in ref. [5], 

\ The major cause of concern is the "hard limiter" effect of the one 

bit processor. With a distributed target such as the sea surface, the dis- 
persed echoes from the various concentric rings "compete" for the quantized 

signal. If there were only 2 reflecting regions, each would only have an 
average receiver output amplitude of % TB (down 6 dB) . If there are 4 signifi- 
cant reflecting rings each would only be % TB in amplitude (down 12 dB) . Thus 
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the desirable property of the 1 bit processor, i.e. that it suppresses clutter 
in an air defense radar, will produce a severe distortion of the impulse response. 

Thus it appears that a "multi-bit" decoder is required. This would 
force a much more complex processor, A study of how many bits are required to 
reproduce the impulse response is given in Part II. While it seems possible to 
get away with a 2-bit I, 2-bit Q system if thresholds are set properly, the com- 
plexity due to pushing the state-of-the-art beyond 1973 technology is disturbing. 
The phase coded system is not recommended at this time. 

3.3 Linear FM Techniques 

The Linear FM or Chirp System is the simplest and by far the most 
widely used form of pulse compression. The primary disadvantage of the ambiguity 
of range and doppler , wherein true range can only be determined when radial 
velocity is known, is not a problem for satellite altimetry. 

The problems of implementation are twofold 

1. Achieving the required bandwidth and dispersion 
with current technology. 

2. Performing integration and digital readout of a 
300-360 MHz signal. 

It is shown in this section that the first problem is not serious 
except that to obtain the desired bandwidth, the only supplier in 1973 is MIT 
Lincoln Laboratory. By the time of actual satellite design, there will be 
many vendors. 

The problem of digitizing a 300 MHz signal with 6-8 bits per word 
is the more severe problem and is the reason for rejecting the relatively simple 
active Chirp configuration shown on Fig. 3. 3,-1, Even if the integration were 
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Expansion/compression portion of a pulse compression 
radar system. 


Figure 3.3-1 




performed with sample-and-hold circuits and an analog integrator. It is not 
believed to be practical. 

Fortunately there are several variations of Linear FM called STRETCH 
or ALCOR that reduce the output bandwidth. These are discussed in Section 3.5. 
There is considerable material in this section on surface wave lines which are 
an essential part of a STRETCH or ALCOR configuration, 

3,3.1 Passive generation of Linear FM signals 

A linear FM waveform may be generated by a passive or an active 
technique. In passive generation, a dispersive delay line is excited with an 
impulse. If the delay line has a bandwidth of 360 MHz, the line output can be 
translated directly by mixing with a transmitter oscillator to the required 
transmitter output frequency. If the delay line bandwidth is less than 360 
MHz, its output frequency may be multiplied to provide the required sweep 
bandwidth, and then translated to the correct carrier frequency. 

The feasibility of a desired delay line is a strong function of its 
dispersion bandwidth product. The state-of-the-art of pulse expansion/compression 
devices is shown in Fig. 3.3-2. For the 360 MHz, 2.8 psec requirement (com- 
pression ratio = 1000), it can be seen that the reflective array compressor 
(RAC) technique described below is the technique (see point A, Fig. 

3.3-2, However, since this is a new invention, procurement would be required 
from MIT Lincoln Lab. Because there are presently no established vendors of 
the RAC line, this approach would involve some development risk if obtained 
from industry. Therefore, an equipment configuration using a lower bandwidth 
delay line followed by frequency multiplication is indicated for the "baseline 
design". A brief description of RAC and other t}^e lines is contained in the 
following paragraphs, 

See section 4,1 for details. 
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DISPERSION - MICROSECONDS 











3.3.2 


Reflective array compressor 


The reflective array compressor (RAC) is a dispersive delay line 
which can be used to provide very high pulse compression ratios at large signal 
bandwidths. This device was originally developed at MIT's Lincoln Laboratory. 

The potential capability of this device is shown in Fig. 3.3-2. The. present 
MIT line has a bandwidth of 50 MHz and a dispersion of 60 |j,s, giving it a com- 
pression ratio of 3000:1. Newer developments are being conducted at MIT and 
elsewhere on shorter lines having up to 500 MHz bandwidth. (Section 3.3, 2.1) 

The technique used in the RAC is an extension of the IMCON dispersive 
delay technique developed at Andersen Laboratories. The basic difference is 
that the RAC uses surface waves instead of bulk waves in the acoustic medium. 

The difference relieves the RAC from the bandwidth limitation of the IMCON 
which is dictated by the thickness of the material, since the desired acoustic 
waves propagate on the surface. 

The RAC represents a significant breakthrough in surface wave dis- 
persive delay lines. The reason for this lies in the fact that the electro- 
acoustic transducers are extremely simple, and are not involved in the dispersive 
properties of the device. The dispersive characteristics are provided by a 
"herringbone" grating etched into the surface of the medium. Previous surface 
wave dispersive lines have transducers which are very large (acoustically) and 
which determine the dispersive characteristics. For this reason, the amount of 
dispersion achievable in a single line has been less than 50 ps, and more 
typically 10 p,s. (A 240 ps line is presently under development). As can be 
seen from Fig. 3,3-2 the RAC is predicted to be capable of dispersions up to 
300 ps. 
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The BtAC geometry is compared with the present IMCON bulk wave 
technique and the conventional surface wave technique in Fig, 3.3-3, Note 
the similarity of the RAC to IMCON, and also the simplicity of the RAC trans- 
ducers compared to the conventional surface wave line. The RAC also requires 
a shorter length of material for the same dispersive delay than the conventional 
design. 

Industry engineers, are following the RAC development closely and 
some already have suggested improvements on the MIT design to overcome some 
of the potential limitations. For lines having the dispersion characteristics 
required by modern radars, these limitations are mainly: (1) acoustic loss, 

(2) temperature sensitivity, (3) spurious responses, and (4) dimensional toler- 
ances. For this radar, the last item is the limiting factor because the wide 
bandwidth forces the line to operate at very short acoustic wavelengths. 


3.3. 2.1 Status of reflective array compressor RAC PC lines 

MIT Lincoln Lab has recently completed 10 RAC surface wave lines 


with the following results: 
Dispersion: 
Bandwidth : 
Weighting: 

Output Pulse: 
Sidelobes : 

Insertion Loss: 

Temperature: 

Weight: 


10 microseconds 

512 MHz max 

Hamming function 

3,5 nsec/ % power 

2 at -25 dB 
others below 30 dB 

55 dB without matching 
45-50 dB with matching 

45° C oven 

Line + matching and shielding 
(without circulators) 

= 0.91 kg 



IMCON 














Temperature: 


+ C for good operation 
Stability: j- .01^ C for exact ranging 

While these lines do not exactly meet the requirements of a satellite 
altimeter, they are close enough, and have the advantage of small weight and size 
as compared to the "multiplier configuration", MIT is willing to supply these 
lines to NASA with appropriate financial support. 

It is felt that this type of configuration would be most suitable for 
a 1976-1980 satellite where weight and power would be a premium. Industrial 
companies (Hughes, etc.) should be able to supply sample lines within 18 months. 
These lines are appropriate for either full analog pulse compression on STRETCH 
techniques discussed in Section 3.5. 

3.3,3 Alternate passive FM generators 

Several alternate approaches to passive FM generation may be con- 
sidered if the bandwidth is reduced. These include, in addition to the RAC: 
IMCON dispersive delay lines, perpendicular diffraction gratings, conventional 
surface wave lines. Each of these is discussed below, 

3 . 3 . 3 . 1 IMCON dispersive delay lines 

Andersen Laboratories, a major dispersive delay line manufacturer 
has delivered "IMCON" delay lines having up to 10 MHz bandwidth, centered 
around a 20 MHz carrier frequency with dispersions of up to 250 microseconds. 

This line utilizes bulk wave propagation in steel. Construction of a line 
of 2.8 microseconds length (see point B, Fig. 3.3-2) should present no design 
problems. The "IMCON" t3q>e line has a linearity of about .01 percent of total 
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phase change. Its cost would be about $10,000 each, with some reduction for 
several units. Thermal control is required for this line to prevent delay 
changes with temperature from affecting system performance. Heater power could 
be held to a few watts by close attention to design of an oven containing the 
line, as well as by controlling spacecraft thermal environment. The line could 
be packaged in about 10 by 10 by 5 cm including heater and oven. Driven with 
an impulse of 1 watt peak, the line would produce an expanded output of -30 dBm. 

The bandwidth limitation of IMCON devices arises from two sources. 
First, the acoustic signal loss in steel increases strongly with carrier 
frequency. Second, spurious propagation modes can occur, if the line thickness 
is greater than one-half of an acoustic wavelength. Since the speed of sound 
in steel is approximately *318 cm per .microsecond, a half-wave length at 
20 MHz is .0079 cm. This thickness of steel is about the minimum which 

can be obtained. The input/output transducers mounted on the edge of the line 
must be bonded very carefully for reliable operation. 

For this radar application, an IMCON line with 10 MHz bandwidth 
would have a time-bandwidth product of only 28. Gating and limiting of this 
waveform would produce considerable distortion of the spectrum which may affect 
measurement accuracy. While TSC has not examined the effect of this distortion 
in detail, it is advisable to avoid it by selecting a line with a larger 
bandwidth. 

3 , 3 . 3 , 2 Perpendicular diffraction gratings (PPDL ) 

A line utilizing this technique is sketched in Fig, 3.3-4. Linear 
dispersion is attained by correct spacing of the transducer fingers. Several 
lines of this type have been in production. Their bandwidth and dispersion 
limits are shown in Fig. 3,3. -2. 
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A line having 3 microseconds delay and 60 MHz bandwidths has been 
built at a center frequency of 120 MHz, and this design could be readily modified 
to the 2,8 microseconds required by this radar. The time-bandwidth product 
of approximately 180 is high enough to eliminate the distortion problems of 
the IMCON, and the bandwidth is low enough to be feasible in practice. The 
perpendicular diffraction grating technique is a recommended candidate. The 
line would be fabricated using fused quartz, and would have dimensions approximately 
2.5 X 3.8 X 1.3 cm. An oven would be required for temperature stabilization. 

3. 3, 3. 3 Surface wave lines 


Several t3q>es of surface wave lines (including the RAC technique) 
are also feasible for signal bandwidths of 60 MHz and dispersions of 2.8 
microseconds (see Fig. 3, 3. -2). While the RAC t}rpe‘ is preferred in terms of 
performance, the conventional designs may be more readily available. However, 
difficulties may be encountered in meeting linearity requirements with the 
conventional designs. 

3.3,4 Status of other surface wave lines 

A status report on surface wave lines for wide bandwidth pulse compression 
systems is given as a result of a visit to Hughes Aircraft, Fullerton, California . 
The primary system that Hughes has built is illustrated on Fig, 3.3-5 and has 
the following characteristics: [4] 

100 MHz bandwidth 
10 microsecond dispersion 
300 MHz center frequency 
30-40 dB insertion loss 
60 dB dynamic range 

Dr. Tom Bristol, Ben Harrington, Hank Gerard (714-871-3232, X4756) 


TB = 1000 
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Figure 3,3-5 






3600 0.1 MIL electrodes 
28 dB sidelobes (single line) 

22 dB sidelobes (double line) 

VSWR 1.3 to 1.5 

Size 7.6 cm by 10.2 cm 

The primary tradeoff in these lines is between the lower losses of 
lithium niobate and the better velocity variation control and lower sidelobes 
of quartz lines. The quartz lines have about 26 dB more insertion loss. 

Lithium Niobate is used in the new RAC lines constructed at MIT. 

With the newer photographic techniques and mesh fabrications, the 
following parameters would be available from Hughes Aircraft Company in the 
late 1973 period. 

200-300 MHz bandwidth 
2-3 microsecond dispersion 
500 MHz center frequency 
40 dB losses (Li. niobate) 

60-70 dB loss (quartz) 

Plus or minus 2^° C yields 1.04 times Hamming 
pulse width, and 0.5 dB loss in S/N 

It can be seen that 1000 to 1 compression ratios are relatively easily obtained. 
However, several companies will have capability for much better performance within 
the next six months to a year. 

There is a procurement out of ECOM, Ft* Monmouth, to build a 250 MHz 
bandwidth line with 40 microsecond dispersion, 30 dB sidelobes, and 50 dB in- 
sertion loss. This is a compression ratio of 10,000 to 1 which is in excess of 
the likely altimeter requirements. Hughes has won this procurement and will 
likely be the first U.S. contractor capable of producing lines with the desired 
characteristics. This, of course, is in addition to the RAC work at MIT. It 
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is likely that Hughes Aircraft will have this capability within one year to 18 
months. Discussions with Raytheon and Autonetics did not yield any additional 
capability. 

These lines are, of course, essential to any Chirp system that might 
be proposed. They would also be used in a STRETCH type system- 


3 , 4 Hybrid Pulse Compression Techniques 

There are several hybrid pulse compression techniques that could be 
considered for altimetry. One possibility is to combine a linear FM ramp with 
a Barker Phase Code. For example, if a bandwidth of 330 MHz was required with 
a 3.3 microsecond dispersion, this could be accomplished with eleven phase coded 
segments of 600 nanoseconds duration. Each segment would then contain a 300 
nsec to 3.0 nanosecond Chirp (TB = 100). The transmit waveform would look 
like Figure 3,4-1-a and the decoded received waveform line Fig. 3.4. -lb for a 
point target. The time sidelobes would not be a problem if the Barker Code 
(length 7, 9, 11, 13) is used. 

The receiver block diagram is shown on Fig. 3.4.-1-C, The signals 
are mixed to a convenient IF and successively pass through a dispersive line with 
a TB 100, a weighting network to reduce close-in sidelobes, and a tapped delay 
line phase coder matched to the Barker Code. 

The advantage of this technique is that the relatively simple pulse 
compression line (TB = 100) is easily made with surface wave techniques. While 
the tapped delay line can also be constructed with surface wave devices at about 
1 GHz center frequency, the tolerances are quite tight. At this point, it is 
felt that the technology will advance within the next two years to make this 
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Figure, 3, 4- la Hybrid of Barker Code and Linear FM 
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Figure 3.4-lc Receiver Block Diagram 


Figures 3,4-1 A Hybrid pulse compression technique. 
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approach slightly more difficult than the "all chirp” system and it would only 
be recommended if surface wave lines of 360 MHz bandwidth with adequate TB were 
not available. 

There are also several hybrids of digital and analog techniques that 
are practical in some circumstances but these generally require processing at 
the full bandwidth of the signal and hence with about 360 MHz bandwidth they are 
not generally attractive. The primary technique that might be applicable is a 
combination of STRETCH and digital pulse compression. If, for example, a 36:1 
stretch were used with a 3,6 microsecond pulse, the received signal for about 
100 nsec of echo would be available over a 3,6 microsecond period and the band- 
width would be 10 MHz. A. digital pulse compression system could be implemented 
using FFT or similar techniques. There is some advantage in that flexibility 
is achieved but at too high a price in hardware complexity, 

3,4,1 Cross correlator 

There is another possible use of the binary phased coded waveform in 
a tracking mode. It comes under the names of "cross-correlator”, "delay lock 
discriminator” and others , A binary phase coded waveform with a known code and 
starting point is transmitted. The transmit code is stored and the code gener- 
ator is started just before the e xpected target echoes (an "early gate”). The 
code is applied to the receiver local oscillator. The mixer output is then a 
decoded pulse the target echo and the delayed code are in coincidence. A 
second delayed code (one bit additional delay) is applied to a second local 
oscillator and mixer for a "late gate”. The difference in the "DC component” 
between the "early" and late gates is the time delay error signal and is used 
as a vernier on the "expected” time delay (altitude) the error signal looks 
like Fig. 3,4-2, 
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Fig. 3.4-3 SWITCHED L.O. LINEAR DECODER AND RANGE TRACKER, RANGE 
DETERMINED FROM CONTENTS OF DELAYED CODE GENERATOR. 
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The processing may be performed at IF or with bipolar video with the 
latter configuration shown on Fig. 3,4, -3. An accurate measure of time delay 
is available on a point target and compensation could be devised for the sea 
echo. The advantages are in the flexibility to change the code and hence the 
resolution and the absence of any dispersive line requirement. 

This technique might be recommended if altimetry were the only goal. 
However, we do not get the full impulse response, and do not see an accurate 
method of wave-height estimation. It will not be considered further. 

3.5 STRETCH-ALCOR Techniques 

A strong candidate for the pulse compression system is the use of 
the STRETCH or ALCOR technique in order to reduce the A/D conversion require- 
ment. In this approach the received signals are partially de- ramped by a linear 
FM de-ramp function which can be generated either by active or passive means. 

The difference frequency signals are now LFM signals with a reduced bandwidth. 

To compress these signals a dispersive delay line is required with a bandwidth 
and delay given by 

TD = T' 

In these equations, 

B = transmitted chirp bandwidth 

SR = desired STRETCH ratio 

N = number of range bins to be stretched 
Rd 

T = transmitted chirp pulse duration 
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The results of the STRETCH operation are compressed pulses whose compressed 

1 SR 

pulse lengths have been increased from ~ to After the STRETCH process, 

envelope detected outputs may now be sampled by an A/D converter operating on 
the reduced bandwidth pulses. This represents one technique for avoiding the 
use of high speed converters operating at the full bandwidth B, 

3.5,1 M ethods of implementation 

There are several ways in which STRETCH or ALCOR can be implemented, 

A general block diagram is shown on Fig, 3,5-1, snd the variations include 

1) Passive Generation (dispersive line) with passive 
generation of a ramp at another slope 

2) Passive Generation with the same slope on receive 
with a filter bank for the "range gates" 

3) Active Generation (swept oscillators) on transmit 
with a passive dispersive line for compression 

4) Active Generation with an active swept oscillator 
on receive 

All of these methods have been tried and the choice depends on the 
flexibility desired and the complexity allowed. For example, the use of passive 
dispersive lines on receive, 2), does not allow a variable slope to examine a 
variable range window. In the active generation technique, more complexity is 
required if the transmit and receive slopes are different. 

One example, that demonstrates the feasibility of building a system 
with the parameters of interest, is the MIT ALCOR system shown on Fig. 3.5-2, 

A transmit waveform of almost 500 MHz bandwidth and 10 microsecond dispersion 
is achieved by multiplying the output of a dispersive line by a factor of 42 
from approximately an 11.8 MHz ramp to a 493 MHz ramp. This ramp is mixed with 
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Figure 3.5-2 Elements of the wideband ramp generator implemented for ALCOR, 
Project Report RDT-13, Lincoln Lab MIT, 20 July 1967. 






a local oscillator for transmission and another local oscillator for a "refer- 
ence ramp", A resolution of 0.5 meters was achieved and sidelobes were expected 
to be in the 30 to 35 dB range. 


3.5.2 


Dispersive line and A/D requirements for STRETCH 


The design equations for STRETCH were used to compute the required 
dispersive line bandwidth (BW) and time delay (TD) required as functions of 
the STRETCH ratio (SR). Results were obtained for a range window correspond- 
ing to 60 contiguous range cells. For B = 360 MHz and T = 2.8 p,sec, the curves 
of Fig. 3.5-3 were obtained. Also, the required A/D conversion rate f^ vs 
STRETCH ratio was obtained, and is shown in Fig. 3,5-4, The parameter K 
represents the oversampling ratio of the compressed and stretched output pulses. 

K = 1 represents taking one sample per output pulse width. Since the compressed 

1 SR 

pulse width before STRETCH is -r , the pulse width after STRETCH is — , For K 

ij B 

samples per output pulse, we get 


If it is desired to keep f^ below 1 MHz (to simplify A/D requirements) , then 
one should use a STRETCH ratio of about 400-1000, 

One might assume that for SR = , the samples could be taken as 

slowly as desired. Note, however, that in this case of complete deramping, it 
is required that range cells be separated by filters of bandwidth 1/T. If 
these filters are to be formed digitally, the deramped spectrum must first be 
A/D converted. Since for the baseline parameters this spectrum has a total 
bandwidth of 11,5 MHz (30 range cells, separated in frequency by .385 MHz for 
T = 2,8 [xsec) , a sampling rate of f^ = 11,5 MHz for the A/B would be required 
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for this case. 


If analog filters are used, the A/D conversion takes place after 
filtering and (with multiplexing) the sampling rate can be less than 1 MHz. 

3,5.3 R ecommended technique s 

An examination has been made of the following cases: 

1. Dispersive line pulse compression with SR = 1 
(Both passive and active pulse generation). 

2. STRETCH pulse compression for SR > 1, followed 
by A/D conversion. (Both active and passive 
generation of the transmit and the receive 
de ramp s igna Is), 

3. Full deramp followed by a filter band to separate 
range returns. (Both analog and digital filter- 
ing are being considered) , 

The first case is just linear FM without STRETCH and has been 
described in Section 3,3. This method was rejected because of the A/D re- 
quirements in digitizing a 300 MHz signal. The choice between the second and 
third case is made on the basis of simplicity. STRETCH pulse compression for 
SR > 1 requires the generation of 3 different linear FM waveforms (one for 
transmit and two for receive deramp signals, see Fig. 3,5-1) whereas full de- 
ramp can be achieved with a single dispersive network. The accuracy and 
tolerance problems associated with matching three linear FM waveforms are an 
added incentive for rejecting STRETCH pulse compression for SR > 1 (Case 2), 

For the full deramp case the analog filter bank is preferable to digital 
filtering, because of the reduction in bandwidth requirements on the A/D converter 
ns described in Section 3.5.2, 


- 71 - 



In summary then, the recommended pulse compression technique is a 
full deratnp method followed by an analog filter bank to separate range returns. 
Considerations concerning the implementation of such a scheme are presented in 
the next section. 
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4.0 


IMPLEMENTATION 


In this section recommendations are made concerning the implementation 
of the selected pulse compression technique. Included are the method of FM genera- 
tion and range processing. An analysis showing the accuracy with which the ramp 
and derarap functions must be generated is also given. 

4 , 1 Recommended Approach for FM Generation 

The state-of-the-art of pulse expansion/compression devices was des- 
cribed in Section 3.3. For the 360 MHz, 2,8 )j,sec requirement (compression ratio 
= 1000) , it was shown that the reflective array compressor (RAC) technique is 
the preferred method (see point A, Fig. 3.3-2.) However, since this is a new 
invention, procurement would be required from MIT Lincoln Lab. Because there 
are presently no established vendors of the RAC line, this approach would involve 
some development risk if obtained from industry. Therefore, an equipment con- 
figuration using a lower bandwidth delay line followed by frequency multiplication 
is indicated for the "baseline design". 

The approach recommended for FM generation is to use a passive dis- 
persive line with a bandwidth of 60 MHz. The vendor of the line would be free 
to use any of the techniques discussed, depending up6n his experience and design 
capability. While TSC would prefer to see the implementation using RAC, PPDL, and 
conventional surface waves in that order, all approaches should be capable of 
meeting the specifications. 

The delay line output is translated to a higher frequency by mixing 
with a 2,33 GHz carrier and selecting the upper sideband output. This operation 
reduces the fractional bandwidth of the FM signal and prevents spectral overlap 
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in the ensuing X6 multiplication process. Standard multiplier designs as well 
as several commercially produced multiplier assemblies are available. The out- 
put of this multiplier chain will be a maximum for a 2.8 microsecond long linear 
FM ramp of 360 MHz bandwidth centered at 14.7 GHz. An RF switch directs the ramp 
to either the transmitter or receiver chain in response to a control signal from 
the central timing unit. The dimensions of the line would be approximately 2.5 x 
3.8 X 1.3 cm plus an oven. 

4 . 2 Range Processing 

In the full-deramp system, a contiguous comb filter set is required for 
processing of the gated and shaped return signal. Total signal bandwidth of the 
comb set is 11.5 MHz and at least 30 range bins are required within this band; 
maximum bandwidth per filter, therefore, is 385 kHz. Center frequency is open 
and may be chosen to give optimum filter performance. 


Four possible filtering methods are available: 

1. Passive LC filters 

2. Active filters 

3. Discrete crystal filters 

4. Monolithic crystal filters 

Features of each of these approaches are listed in Table 4.2-1. 

The requirements placed upon the filter set are rather minimal. The 
quality factor, or Q, which is the primary determinator in filter feasibility, is 
only 50. Input signal level can be made sufficiently high that ultra- low filter 
losses are unnecessary. Furthermore, since only 30 range bins are required, size 
is not a major factor. 
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TABLE I, Features of Filter Approaches 



PASSIVE LC 
FILTERS 

ACTIVE 

FILTERS 

DISCRETE 

CRYSTAL 

FILTERS 

MONOLITHIC 

CRYSTAL 

FILTERS 

Power Consumption including 
Input Amplifier 

LOW 

HIGH 

LOW 

LOW 

Siz:e 

Large 

Larger 

Medium 

Small 

Reliability 

Good 

Many passive 
components 

Fair 

Many active 
and passive 
components 

Good 

Few components 
Proven crystal 
reliability 

Probably Good 
Few components 
Monolithic 
Construction but 
unproven 

Calibration Effort 

Large 

Many Adjustments 

Large 

Many Adjustments 

None 

No Adjustments 

Small 

Only Frequency 
Shift requires 
Adjustment 

Parts Cost 

Low 

Medium 

Medium 

Highest ' 

Possible 

Modifications of 
Monolithic Masks 

Design Risk 

Low 

Medium 

Low 

High* 


Both risk and cost of monolithic filter may be smaller in 1980 time-frame 



To work at the rather wide bandwidths (about 2%) required in this system , 
crystal filters, both monolithic and discrete, would require paralleling. Since 
crystal filters at the required I.F. are limited to approximately 100 kHz band- 
width, four filters would be paralleled for each range bin. This would be possible, 
but the advantages of crystal filters -- size, weight, and reliability -- would not 
be realized . 

Both active and passive filtration would be feasible at these frequencies. 
Indeed, to a first approximation, active and passive filters have the same A-C 
characteristics. Active filters, however, would consume power. If each active 
filter were to need only 100 milliwatts of D-C power, then the comb filter set 
would require about 3 watts. While this figure is entirely reasonable for a ground 
or airborne system, it is not at all desirable in a spacecraft. 

The contiguous comb filter set baseline design, therefore, is chosen as 
passive b-C filtration. A major disadvantage of this scheme is that adjustments 
would be required. Calibration time for the airborne system would be relatively 
large, and the spacecraft system would require further effort in either selecting 
fixed components to replace adjustable ones or in fixing the adjustments inmovably. 
Filter reliability would be excellent, however. Only low-failure-rate passive 
components would be required, none of which would be electrically stressed to any 
significant degree. 

The contractor would have the option of either building these filters or 
buying them from an outside company. Many manufacturers make L-C comb sets as a 
standard product line. Alternately, many aids exist for internal filter design. 

Following the filter bank, detectors and a means of sampling the detector 
outputs will be required. Simple diode detectors may be used on the filter outputs; 
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analog multiplexers will then follow these detectors. Inaccuracy due to droop on 
the detector outputs may be circumvented by using FET multiplexers with very large 
OFF resistance and following the multiplexer output with a fast sample-and-hold 
amplifier. This scheme allows extended sampling of the filter outputs during as 
much as the full pulse repetition interval, and therefore permits very light, low- 
cost analog-to-digital conversion for data transmission. All devices used in the 
detector are light-weight and low power units. Most are made as integrated circuits. 
Indeed, many such devices are already NASA-qualified and/or MIL-qualified . In 
addition to the units already described, some timing and logic will be required 
for clocking the multiplexers and the A/D conversion chain. 


4 . 3 Ramp /De-Ramp Generation Requirements 


This section will present an analysis of the accuracy with which the ramp 
and de-ramp functions must be generated. It must be kept in mind that since pulse- 
to-pulse processing is not employed, the accuracy requirement need only be imposed 
on a single pulse basis. After a ramp is generated (either actively or passively) 
and transmitted, a de-rarap signal will be required within a time interval less than 
1 millisecond which matches the ramp function within some accuracy. The following 
analysis will determine this required degree of matching of the two functions and 
the linearity requirement which both must achieve. 

Consider the transmit pulse which ideally would be a linear FM signal 
of duration T seconds and bandwidth B Hz. We will model the actual ramp as a 
signal whose frequency versus time is given by 

f^ = f^t + sin U)t 


where 

f^ = slope of ramp 

Afj^ = peak frequency error (deviation from ideal ramp) 
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(jj = modulating frequency of error 
ul)T 

Note that 't" = number of cycles of sinusoidal frequency error across the 
2tt 

transmit pulse. 

Next consider the de-ramp signal. This signal will be initiated at 
the expected time of the first return. Since the last target return of 
interest will be received a maximum of At seconds after the first, ideally 
the de-ramp signal duration should be T + At seconds. However, since 
T = 2.8 ysec and At~140 nsec, the de-ramp signal could also be T seconds 
long with little loss of received energy. The de-ramp signal can be 
written 

f 2 ~ ^2 + Af (t - At) + Af^sin o)(t - At) 

In this equation we have allowed the de-ramp to differ from the received 
signal in four ways: 

1. it can have a different slope, to allow for partial de-ramp (STRETCH), 

2. it can have a slope error Af, which could arise if the ramp and 

de-ramp were generated by separate techniques or separate circuits, 

3. it can have a different amplitude of frequency error, ^f 2 » 

4. it can have a delay error At which will vary from zero to about 140 

nanoseconds, for the first and last target returns respectively. 

After mixing the received signal with the de-ramp signal and retaining 

only difference frequency terms, we obtain the following signal 

• • • • * 
f^ - f2 = (f^-f2)t + ^2^^ “ 

+ [Afj^ sin (jjt - Af 2 sin o)(t-At)] 

These five terms will be analyzed separately to determine their effects on 
the range response of the pulse compression system. 
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First Term 


This term is a residual LFM term which will be required whenever 
STRETCH is used. The slope is assumed to be precisely matched to the 

STRETCH dispersive line. The effect of a mismatch is included in the 
third term and will be analyzed there. 

Second Term 

The term represents the desired frequency shift due to the range 

difference between the first target return and one displaced from it by a 
distance corresponding to At seconds. In the case of full de-ramp, for 
example, this term represents the conversion of range displacement to 
frequency displacement. Adjacent range cells are thus detected separately 
in the outputs of adjacent filters. 

Third Term 

If the slopes of the ramp and de-ramp functions are in error, the 
signals after de-ramp will contain an error frequency term which is linear 
with time. This gives rise to quadratic phase effects which degrade range 
resolution and range sidelobes. Exact computer simulations of the system 
impulse response show that one can allow about tt/ 2 of quadratic phase to 
build up from center to end of the pulse duration without significant 
degradation of the range resolution or sidelobe levels. 

Slope Accuracy of Generating or f ^ 

The accuracy is determined on the basis of allowing tt/2 quadratic phase 
to accumulate (center- to-end) due to an error Af. Consider a signal 

* • f¥i m 

f =* (f + A£)t 7 ^ ^ ^ 2^ 
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Then the quadratic phase error due to Af is given by 


T/2 . 

A4> = 27T ^Af t dt 

Jo 


Af T 


Setting 


A(j) = tt/ 2, we find 



Note that f = 

T 

Expressed as a percent error in slope f, we find from the above 


M = _2 

^ BT 

For B = 360 MHz. T = 2.8 ysec, 

^ 0 . 2 % 

f 

Total excursion Af is given by 

Af = Af T = I - .714 MHz 

Hence, if we wish to sweep from zero to 360 MHz, the end point frequency 
can be 360 + 0.714 MHz. 


- 80 - 



Fourth Term 


The term Af At represents a constant frequency shift which will give 
rise to a range measurement error. In order to make its effects negligible^ 
we can require that the frequency shift correspond to no more than one- 
tenth of a range resolution cell. Since one range cell translates to a 
frequency shift of — Hz in the case of complete de-ramp (i.e. , adjacent 
filters are spaced by Hz), we will impose the constraint 


Af At < 


- lOT 


Note, however, that the quadratic phase constraint has already led to 
the relationship 



Hence using this allowable value of Af in the above gives rise to 
^ Por our baseline design, At = lAO nanoseconds and T = 2.8 ^seconds. 

Hence , we conclude that if the quadratic phase constraint is met, the fre- 
quency shift due to the fourth term will not degrade the achievable range 
measurement accuracy significantly. 


Fifth Term 

The fifth term, due to an assumed sinusoidal deviation of the ramp or 
de-ramp function from ideal linear FM waveforms, can be rewritten in the 
form 


f = Af sin (tot +0) 
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Where 


Af = /(Afj^ - Af 2 cos wAt)^ + (Af 2 sin toAt)^ 


6 = tan 


Af 2 sin u)At 
Af^ - Af 2 cos o)At 


Examination of the above will show that if Af^ ~ Af 2 » then Af will range 
between zero and 2Af^. With little loss in generality of results, we can 
ignore the constant phase angle 6 and consider the worst case example of 
a frequency error 


f = Af sin wt 


Sinusoidal Frequency Error 

The phase error associated with this frequency error is given by 


A4> = 2tt /* Af sin o) t dt = cos wt 


The peak value of this phase error is 


. . 2irAf Af 
A<J) = = — 

P w f 


f = modulation frequency, i.e» rr- 
m ^ * 2 it 


For a system utilizing full de-raraping, we integrate the de-ramped 
signal for a time T seconds, to form a filter of -4 dB bandwidth This 

filter has a sin x/x shape with its first null at Adjacent filters are 

1 1 

spaced ~ Hz, in order to compress pulses separated by At = — (see sketch 
^ B 

below) : 
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For At 




If after de-ramping we get a CW pulse of length T which has a small resi- 
dual sinusoidal phase shift of peak value A(|> due to a sinusoidal frequency 

P 

error in the de-ramping signal, then paired echos in the filter output 

A(j) 

response will be created of value — relative to peak response and displaced 

from it by + f the modulating frequency of the error. For one cycle across 
“■ m 

1 N 

the de-ramped pulse, f = “. For N cycles, f “ If we wish to keep this 

mi ml 

sidelobe response (which is equivalent to a range sidelobe) less than 23 dB, 
for example, we need 


A4> 

20 loe — 2- 


-23 


so 

A(f) = 0.14 
P 


Hence 


Af * 0.14 f 

m 
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For f = — (probably the lowest, or worst case value), 
m T 


Af = 


0.14 


0.14 

2,8 X 10 


-6 


= 50 kHz 


Since in the worst case, Af = 2Af^^ we obtain Af^ = 25 kHz, This is the 
allowable peak sinusoidal frequency deviation from a straight line for the 
case of one cycle variation across the pulse. For higher frequency variation, 
the allowable peak error increases. For 10 cycles across the pulse, 

Af = 500 kHz, and Af^ = 250 kHz. 


4,3,1 Suimnary 

The following summarizes the conclusions reached regarding the required 
accuracy of generating the ramp and de-rarap functions. Those conclusions 
apply to either active or passive generation techniques. 

1. The ramp and de-ramp functions must have slopes which meet the 
following linearity requirements: 
a. Full De-Ramp 


fl - f2 1 Af 


b . STRETCH 




where f^. = design slope for selected dispersive 
^ STRETCH line 


In both cases, the tolerances on Af are given by Af = —j. For the 
selected design (full de-ramp) this leads to a percent linearity 
requirement of 


K 

f 


< 0 . 2 % 


- 84 - 



2 . 


If (1) above Is met, then the range measurement error due to 
♦ 

the Af At term will be negligible for the selected baseline system. 

3. To keep slowly varying frequency deviations from degrading sidelobes, 
the peak allowable frequency deviations are: 

a. Af^ < 25 kHz (one cycle of variation across pulse) 

b. Af^ < 250 kHz (ten cycles of variation across pulse) 

When “ < 1, these very slowly varying frequency deviations begin to look like 
2tt 

slope errors and are essentially accounted for by the specification on Af in 
(1) above. 
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PART II 


SUPPORTING STUDIES 
1 . 0 ^INTRODUCTION AND SUMMARY ' 

In Part I, the design of a high resolution satellite altimeter was dis- 
cussed. During the design effort, several supporting studies were conducted. 

The results of those studies are reported here. 

Section 2 discusses the investigations of improved altitude tracking 
algorithms. The form of the maximum likelihood estimator (MLE) for simultaneously 
estimating epoch, signal- to-noise ratio, and wave height is derived, and its per- 
formance is computed. In addition, modifications to the standard split-gate 
tracker are analyzed, and it is shown that the split-gate tracker performance 
can be improved by increasing the width of the late gate and lowering the position 
of the early gate on the leading edge of the return. 

Section 3 considers the performance of a digital binary phase code 
receiver. It is shown that a multi-bit digital processor of at least 2 bits plus 
sign is required for adequate performance. This result is due to the distributed 
nature of the target (sea surface). The results show that for a compression 
ratio of 1000 and an uncompressed signal-to-noise ratio of -20 dB, the compressed 
output signal-to-noise ratio is only 4.1 dB. This is a 5.9 dB loss relative to 
a linear system. A system with 2 bits plus sign gives an output signal-to-noise 
ratio of 7.9 dB (i.e. about 2 dB loss relative to a linear system). 
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Section 4 makes a comparison of Barrick’s model of the sea surface return 
with several experimentally derived values. It is found that Barrick’s model agrees 
reasonably well with the data trend versus wind velocity. However, it is found 
that his model predicts values for the cross-section per unit area which are some- 
what higher than the measured values, and a reduction of Barrick's model by about 
6 dB gives a more conservative model. 


Section 5 describes the recommended changes in the system parameters for 
an aircraft mounted altimeter. It is necessary to modify the selected set of 
system parameters given in Part I, since that design is based on orbital parameters 
whereas the breadboard unit will actually be tested using an aircraft. 
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2.0 


IMPROVED ALTITUDE TRACKING ALGORITHMS 


2. 1 I ntroduction and Sunanary 

The general problem in altimetry is that of estimating the time of 
occurrence, epoch, of the leading edge of radar reflection from an extended 
noise-like target., (i.e., the sea surface). Here, a variety of techniques 
for estimating this epoch are developed and compared with the performance of 
a standard split-gate tracker. Both the optimum, maximum likelihood estimator, 
and sub-optimum, "modified split-gate" tracker, techniques are shown to provide 
considerable improvement in performance when compared with the conventional 
split-gate tracker. 

In Section 2.2, the forms of the maximum likelihood estimator (MLE) and 
of the minimum mean square error estimator (MMSE) are derived. The as 5 mipotic 
variances of the MLE for simultaneously estimating epoch, wave height and signal- 
to-noise ratio are computed. 

For the special case in which the wave height is assumed known, the 
MLE and MMSE arc shown to be very similar to a split-gate tracker. In both cases, 
the width of the late gate is much larger than the early gate width. It is also 
shown that at high signal- to-noise ratio the MLE places the center of the early 
gate well below the half power point on the leading edge of the return. 

These results suggested that a conventional split-gate tracker which has 
equal early and late gate widths and which centers the early gate on the leading 
edge of the return signal may not be optimum. By appropriately designing the 
gates and their position ,the performance of a split-gate tracker may be improved. 
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In Section 2.3 the performance of a split-gate tracker is derived, and 
the effect of changing the gate widths and the placement of the early gate is 
examined. For the parameter values considered, it is shown that the standard 
deviation of the tracking error can be reduced by a factor of 0.74 by increas- 
ing the width of the late gate. If, in addition to increasing the width of the 
late gate, the early gate is centered on the quarter power point of the leading 
edge of the return, the standard deviation of the tracking error is reduced by 
a total factor of 0.5. 

The performance of an "optimized" split-gate tracker is computed. The 
optimized tracker shows significant improvement over the conventional tracker. 

In particular, it does not display the saturation effect at large signal-to- 
noise ratios. However, the optimized tracker must be adaptive in the sense that 
it has to vary the position of the early gate as a function of signal-to-noise ratio 

Bias errors of a split-gate tracker are computed and shown to be a 
function of wave height and signal-to-noise ratio. In general, the tracker 
output must be corrected for these errors. It has not been determined whether 
these errors are more easily compensated in the tracker itself, or by "off-line" 
QQ^j- 0 Ctions on the ground. It should he noted however, that a tracker which 
compensates for signal- to— noise ratio and wave height is, in effect, estimating 
these parameters as well as epoch. Such a tracker may be nearly as complicated 
as a MLE, 


In the initial work on the split-gate tracker, it is assumed that the 
data is sampled at time intervals far enough apart in range so that the sampled 
values are independent (descrete case) , Since this is an unrealistic assumption 
for an analog split-gate tracker, the split-gate tracker equations are re-derived 
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assuming that the sampled values are no longer independent. The continuous 
case (analog) is then approximated by letting the time samples occur at very 
closely spaced intervals. 

In Section 2.4, the accuracy expressions derived for the continuous split- 
gate, the discrete split-gate, and the MLE trackers are evaluated and compared. 

2 . 2 Optimum Processing 

From standard statistical theory, it is known that the Cramer-Rao bound 
gives a lower bound for the variance of unbiased "regular" estimates, and (if 
there are a sufficiently large number of samples v^) maximum- likelihood processing 
approximately achieves this lower bound. As usual, the derivation of a maximum 
likelihood processor for altimetry begins with an appropriate set of assumptions 
and definitions. 

Here, it is assumed that the received signal is a stochastic process 
which is Gaussian, has zero ensemble mean, and is nonstationary: specifically, 
its average power varies with respect to time. It is also assumed that 

(a) The average power of the signal stochastic process can 
be described by a function f(t, a) where a is a multi- 
dimensional parameter. 

(b) The observed data actually consist of a set of samples 
of the received p ower at a discrete set of times t^, 

(c) The objective is to estimate some component of 0^, 
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where 


Of particular interest is the case where 


f(t, a) = a g (t~t^, 


g(t.-t^, 3) is the normalized signal shape (mean power) 
returned at time t when the epoch is t^. 


1 


is a vector of parameters needed to describe 
the return signal shape (e,g. wave height, 
sidelobe level etc.)* 


and 

t is the time of arrival (epoch) of the signal, 

o 

a is the signal-to**receiver-noise level. 

If it be assumed that the received power is measured at a set of time 
instants t. separated by a range resolution cell width of the radar, then the 
observed data consist of a set of samples which have the following statistical 
properties : 


(a) V, are statistically independent 

(b) has probability density functions 


P.(V.) dv = — exp [-v./v.] dv 


where 


V. is the returned power (relative to receiver- 
1 

th 

noise) from the i — signal sample. 


and 


V. = E(V.) = a g(t.-t ) + 1 
1 1 r o 
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Now, as is well known, maximum likelihood processing would consist 


of maximizing the joint probability density function for a set of N samples, 

N 

p(Vj, V 2 . . . .) =Tl (l/'V.) exp (-V./V.) 

i=l 

The set of N resolution cells, i=l,.,,,N, is a specific set of cells-- 
say, i=l corresponds to the cell at some specified range beyond the "first target". 

Maximization of P is the same as minimization of -log^ P, Thus, the 
maximum likelihood estimate of t^ is obtained by minimizing with respect to a, 
t^, and ^ the quantity 


^MLE^^o' 




N 

= S 
i=l 


V, 

1 


ag(ti~t^,J) + 1 


N 

+ S log^ [ag(t.-t^,^) + 1] (2.1 

i=l 


If the minimizing values are given by a, t^, then the maximum like- 
lihood estimate of t is t • Furthermore a and 6 are also maximum likelihood 

00 

estimates of a and _g, respectively. 
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It is also well known that the minimum mean square error is obtained 
by minimizing the function: 


1=1 


( 2.2 


Thus the form of the MLE or MMSE processor can be obtained by differ- 
entiating ^MMSE respect to the appropriate variable, and setting 

the result equal to zero. 


2,2,1 Form of the MLE and MMSE 

As a first step towards obtaining the MLE and MMSE, it will be assumed 

that the signal-to-noise ratio, a, the signal epoch, t^, and the waveheight, 

0 , are unknown. The shape of the function for the returned power will be that of 

[ 1 ] 

an ideal receiver and the range sidelobes will be neglected. Thus, 

g(t,^) = POt) (2.3 


where 


t 


dU 


(2.4 


is the normal probability integral and 


3 
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rms waveheight 


c = velocity of light. 


Now let 


V(t) = ag(t-t^,3)+l. 


Then 


N /V. 
1 


1=1 \ V 


(2.5 


and the MLE for t^,a,3 must satisfy the equations 


^^MLE ^ /^i - V \ av 

at .1 2 I at " ^ 

° " \ V / ^ 


( 2.6 


^^MLE ^ 2 / ^i - V 


a a 




V 


(2.7 




MLE 


^3 


= S 

i 


^i - V \ a^ 


V 


a3 ° 


( 2.8 


Now 


a V / 7 7 

^ exp/-3 

o 2n 


(2.9 
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(2.10 


^ = -^ = ”75^“ 



(t.-t J 


( 2.11 


And the MLE must satisfy: 



V. - ap(e(t.-e^))-l 
(^ap(e(ti-to))+i)^ 



(2.12 


and 



V. - ap(3(t.-t^))-l 
(ap(3(t.-t^))+l)2 


p(3(ti-to))= 0 


(2.13 



V. -ap(e(t,-t^))-l 
(ap(3(t.-to))Pl)' 


(t.-t )exp(- 3?3^(t,“t )^) = 0 
L o \ 1 o ' 

(2.14 


The equations in this form do not give much insight into the form of 
a MLE; however, if a simplifying assumption is made, the nature of the MLE 
becomes more apparent. 


The first assumption, is that the sample rate is high compared to the 
rise time of the leading edge, (i.e., the Nyquist rate or higher). Then, 
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P(3(ti-to))+l/a ^ 0z\t 


exp(-%3^(t^-t^)^)3At 
F(3(ti~t^)) +l/a 


2 

I f exp(~^u ) du == ^>[LE 
3At y P(u)+l/a 3At 


(2.15 


and similarly. 


!l!lVV)_«^ r P(u) du = >4E 

P(3(fci-to))+l/a PAt y P(u)+l/a 3At (2.16 

i 

and 




i 


(t ~t ) exp f-%3^(t -t )^) 
1 o V X o / 

P(3(ti-to))+l/a 


1 ^ u exp (-%u^) 

P^AtJ P(u) + 1/a 


du 





(2.17 


and equations (2,12 and (2.13 become: 


E ''i 4 e = 0 
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r V. (t.-t^) h(3(t.-t^)) 


a 


(2.20 


e^At 

where; 

[P(t)+l/a] 

and 

(P(t)+l/a)^ 

These functions are plotted for several different signal- to-noise 
ratios in Figs. 2.1 and 2.2, As might be expected, at higher signal-to-noise 
ratios, the weighting functions place more emphasis on the leading edge of 
the pulse. A more surprising feature is the relatively high weight placed 
on signal values which are equal to or less than the noise level. 

Heuristically, the MLE estimators may be considered as being implemented 
with three filters. The shape of each filter's response is governed by the 
signal-to-noise ratio and the time scale factor is governed by the wave- 
height parameter, A brute force approach to obtaining the MLE would then 
be to repeatedly filter the data, changing the parameters a and 3 on each 

pass until an epoch, t^ is found for which the equations 2.18 - 2.20 are 
satisfied. 


( 2.21 


( 2.22 
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If equation 2.2 is differentiated, the MMSE is obtained. As it turns 
out, the form of the equations to be solved is identical to that for the 
MLE. However, the constants and filter impulse responses become: 




4iSE = JO>(«)+l/a)P(u) du 


■^MSE ° 


du 


(2.23 






The MMSE weighting functions are plotted in Figure 2.3. 


If the wave height, or equivalently the returned pulse shape is assumed 
known, then the MMSE for epoch becomes equivalent to a split gate tracker. 
That is, since 3 is assumed known,’ then only eq, 2,18 and 2,19 need to be 
solved. The signal-to-noise ratio, a, can be eliminated from these equations 
to yield: 
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IMPULSE RESPONSE 


MMSE Filter Impulse Response 


u “ 3t 

s 

max 
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S V. h (3(t.-t^) - C r V. g (3(t.-t^) - B = 0 
i ^ i 


(2.24 


where 

' 2 
C = JP(u)exp(-^u 

JP^(u)du 
B = + C 

Thus the estimate for epoch is obtained by passing the data through 
an early gate with impulse response h(3t) ; subtracting the output of a late 
gate with impulse response Cg(Pt); subtracting a bias term B, and then 
taking as the epoch, the point at which the resulting output passes through 
zero. 

2,2,2 Variance of the MLE 

The asympotic variance (as the number of pulses averaged increases) 
of the MLE is computed in this section. Since all the trackers require a 
large number of pulses to be integrated for the design goals to be met, it 
is felt that the asympotic analysis is appropriate. 

It is known from the theory of maximum likelihood estimators (see 
Wilks [ 2 ]) that the covariance matrix for the MLE is the inverse of the covari- 
ance matrix of the score. 
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The score is defined is follows: 


be the gradienr operator. Then the score is defined as 


S(Vi. v^. . . = V log P(V^, V^. ... 

Thus the covariance matrix for the MLE is given by 


"mle < 

Lank [3] shows that for the altimetry problem, the covariance matrix of the 
score is given by: 


rr vT 


= ^ (V log V^) (v log V^) 




Recall that 


V. = 

X 


a P ( 3(t^-t^) + 1 


so that 


P (6(t^- t^) 

-a6 P _ (S(t^- t^) 
a(Kto) P (6 <t.- t^) 
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Thus 


Where: 


and where 




P. 

1 


POCt.-t^)) 


( 2.30 


( 2,31 


( 2.32 
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The computation of the above matrix can be simplified if it is assumed 
that the samples are closely spaced in time: Then 


R 


aa 




a^3At 


^ Sit 


or 


R 

aa 


N 

P 

a^3At 


/ 


F(u) 

H(u) 


du 


(2.33 


(2.34 


where 


H(u) = P(u) + 1/a 


(2.35 


Applying similar approximations to all the correlation coefficients yields: 


K 

^aa " . N 

2 P 

a 3 At 


R. 


at 


R. 


h3 


t t 
o o 


‘t 3 
o 


R, 


33 


aAt 


K 

■a3 . N 

2 P 

a3 At 


A 0^0 P 

At 


1 K 
3At 


t 3 • 

o P 


bV 


2 ^ *^8 . N 


aa 


=/[^1 

[h(u)J 


du 


^ato . N ; K.^ = f P(u)PCu)_ du 

-Am. “ -J *> 


H^(u) 


K33 = j uP(u)P(u) du 


H^(u) 


[h(u) 


:] 


t 8 =/“ 

° [h(u)J 

86 “I 


du 


uP(u) ^ du 
H(u) 


(2.36 
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The coefficients can be evaluated directly using: 


P(u) 


^ 1 




./ 


^ 00 


exp(-t^/2)dt 


(2.37 


However to do so^ numerical integration techniques must be used. To simplify the 
calculations, P(u) is approximated by 


P(u) = au + 1/2 


u e (-l/2a, l/2a) 


1 


u e (l/2a, ') 


0 


u e (-“ , - l/2ot) 


(2.37 


The constant, a, is chosen so that the ramp defined above approximates the 
normal integral in a minimum mean square sense. Thus: 


.3227 


(2.39 


Using this approximation yields 


K 


aa 


= 1 \1 - 2 los (1 + &) +( 1 \ + (3aT^^^-l/2) a 

\l+a/ 


max 


(1+a)^ 


K =-jlog (1+a) - a 

^ o L 1+aJ 


K 


a3 


1 f 1 

a 


/ a+4 \ log (1+a) + (a-i-2) 

1 2a / 2 (1+a) 


t t 
o o 


= aa 


1+a 
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yields 


/ N D KD 


So that 

R 

MLE 

In particular, 
0 * 


= R 


-1 


gAt d“^ 

N 



-1 


K 


-1 


1/2 

aa 

1/2 


t t 
o o 


and 



( 2.40 

( 2.41 

( 2.42 

( 2.43 

( 2.44 

( 2.45 

( 2.46 

( 2.47 
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Note that matrix K depends only on the parameters 

U ag T (2.48 

max max 

and a. Here T is the time measured from the half power point to the end of 
the observation interval. 

—1 ^ —1 ^ —1 ^ 

The quanities [K ] , [K ^ > and [K ] are plotted in Figure 

aa L t _ _ 

o o gg 

2.4, 2.5 and 2.6. 


- 108 - 






S/N In dB 
- 110 - 









2,3 


Split-Gate Trackers 


In this section the performance of a split-gate tracker is derived, and 
the effect of changing the gate widths and the placement of the early gate is 
examined. Bias errors are computed and shown to be a function of wave height and 
signal-to-noise ratio. Finally, the model is extended to include correlated 
samples to approximate an analog (continuous) tracker. 

2.3.1 Discrete Case 

2.3. 1.1 Variance of a split-gate tracker 

To keep the derivations relatively simple, a sampled system will be 
assumed, A split-gate tracker estimates the epoch of a signal by balancing the 
output of an early gate with that of a late gate. Specifically, let: 

t^ , i = 0, +1, +2,. , .be the times at which the signal is sampled. 


v(tp 


S 6(t .-t)h, 
11 

1 

C^(ti“t)gi 

i 


be the square law detected output of the receiver 

at time t . , 

1 

be the impulse response of the early gate, 
be the impulse response of the late gate, 


and 


C be a constant (to be determined later) which makes 

the tracker unbiased. 
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Note that the Kronecker delta functions arise from the assumption of a sampled 


system. 


For any time delay, , between the tracker estimated epoch 
and the true epoch, the error voltage of the split-gate tracker is given by: 


e(T^) = 2v(t^-TQ)(h^-cgp 

i 


(2.49 


The tracker feedback loops function so that the error voltage is held at zero. 
Thus the tracker estimated epoch is given by the solution of: 


= 0 


(2.50 


If the loop time constants are long, and the tracking error is 
unbiased, then will be nearly zero. Expanding C^(Tq) in a Taylor series 
about Tq = 0 and substituting into (2.50 yields: 




3 ^ 


T=0 


(2.51 
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Let 


k = E 5£(T) 
ST 


be 


T=0 


the mean tracking error slope. Using (2.51, one has: 


Var 


Var g(Q) 
k' 


(2.52 


The quantities k and Var £(0) depend on the shape of the early 
and late gates, and on the characteristics of the mean power return, V(t) , The 
special case sketched in Fig, 2,7 is now considered. As in Section 2.2.2, it is 
assumed that the ramp is chosen to approximate the normal integral in a mean 



Fig, 2.7 Sketch of the average returned power level, 
and the gate positioning. 
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That is, the following simplifications are made; 

i) Square - non- over lapping gates are assumed, 

ii) The late gate is positioned on the plateau, 

iii) The early gate is positioned on the leading edge, 

iv) The signal return power is modeled as a linear ramp rising 
to a flat plateau. 

v) The receiver noise level is unity. 

If the following definitions are made; 

At = time between samples 
a = signal-to-noise ratio 
T = rise time of the leading edge 

M = T/2At; (2M+ 1)= total number of samples available on the 
leading edge 

N„ = number of samples in the early gate 
£ 

N, = number of samples in the late gate 
1-1 

Vo = average power return at the first sample in the early gate 

N = number of pulses averaged 
P 

It is shown in section 2.3. 1,4 that the tracking variance is 

given by; 
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Note that for a fixed early gate width, N , the variance of the 

hi 

tracking error, is minimized by making as large as possible and by making 
V as small as possible. This is consistent with previous results which 


showed that the MLE used as large a late gate width as possible, and tended 
to push the early gate position as far down the leading edge as possible. 


In general, the optimum width for the early gate depends on the 
ratio Vp /(a/2M); however^ it can be shown that for high signal-to-noise ratio, 

a, the optimum early gate width is ^^,=1. In this case, the tracker variance 

\ 

becomes 


Var Tq = 


( 2Mt) 


N 




( 2.54 


Thus for the "optimized" split- gate tracker, the standard devia- 
tion of the tracking error is inversely proportional to signal-to-noise ratio. 

A consequence of optimizing the tracker gates is that the tracker saturation 
effect no longer occurs. This improved performance does, however, require in- 
creased complexity in the tracker. Note that in order to achieve optimum perfor- 
mance, V , the average power output from the early gate, must be held constant 
independent of the signal-to-noise ratio. It can be shown that this can be done 
only if the normalizing constant is changed as a function of signal-to-noise 
ratio. Thus an "optimized" split-gate tracker would have to be adaptive in that 
it must estimate signal-to-noise ratio, and then adjust the early gate position 
accordingly. 
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2. 3*1. 2 


Conventional Split-Gate Tracker 


In most cases, it is more convenient to design a tracker which 
tracks a power level that is a constant ratio of the plateau level, e.g.^fthe 
half power point or the quarter power point. In this case, the normalizing 
constant does not have to be adaptive at high signal-to-noise ratio. To compare 
such a tracker to the "optimized" split-gate tracker, consider the case: 

O- 

N_ = 2, M = 1, and V = 1+k a (i.e, the tracker is adjusted to track the k 
power point (0<k<l) ). For this case, the tracker variance is given by 

2 

Var (1/32 + %k^(l+l/Y) + k/a(l+l/Y) + (%a"^) (1 + 1/Y) ) 

P 

(2.56 

In the above expression, Y is the ratio of the late gate width to the early 
gate width (i.e. Y = N /N ) , The effect of increasing the late gate width 
(increasing y) ^nd of changing the tracking point, k, is shown in Fig. 2,8. 

The case k = %, Y = 1, corresponds to the conventional split-gate tracker. 
Relative to this case, increasing the late gate width, and lowering the position 
of the early gate both improve tracking performance. As can be seen, the 
greatest effect seems to be due to lowering k (the early gate position) . 

For comparison the "optimized" split-gate tracker performance 
is also shown. As can be seen, considerable performance improvement is possible. 


(1 

Note tha^ Vg i.s the average power at the midpoint of the early gate. In 
general • 
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2.3. 1.3 


Split-Gate Tracker Bias 


A common practice for reducing tracker bias error entails subtract- 
ing an estimate of the noise, from the return waveform, i.e. clamping the 
return waveform to In this case, the normalizing constant (Section 2. 3. 1.4) 

is given by; 


N V 

r = ^ E 

a 1-T| +a 

L 'o 


where V_ is the average power at the midpoint of the early gate. If the 
E 

tracker is adjusted to track the k power point, then 


Vg = + ka 


and 


C 


a 


!e 

\ + a 


(2.57 


If the constant is set for the wrong signal-to-noise ratio, say a^ , then 

there exists a bias t such that the mean error signal is zero at T . Thus, 

o o’ 


T 

E(e(T^)) = 0 = + a(k+ ^)) - Cg Nj_a-T^ + a) 

o 


(2.5S 
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and the bias normalized to waveheight is given by: 


T i-'n 

(I-P) + (k-p) (2.59 

where 

1-T] + ka 
o o 

p = (2.60 

Note that while the bias is a function of waveheight and signal- 
to-noise ratio for any k 1, it can be made arbitrarily small by accurately 
estimating the noise level. That is, as the bias approaches zero since 

p -»k. 

2,3. 1.4 Derivations 

The tracker error signal is given by: 

ecy = i: j^v(t . -Tg) (h. -c g^j (2.6X 

The tracker will be unbiased if 
E G(0) = 0 


Thus the normalizing constant is given by: 

S V(t.) h. 

C = ■ — i ^ (2.62 

2v(tp g. 
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where V(t^) = E V(t^) is the mean receiver power output at time t^. 


The mean error slope is given by 


dE5(r: 


- V(t.)(h.-c g^) 


(2.63 


And the variance of the tracking error signal is given by: 


Var e(0) = 2 Var[v(tJ (h^^-C g^' 


(2.64 


where the samples are assumed to be uncorrelated, (i.e. E(V^-V^) (Vj-Vj)=0, i^j) 


For Gaussian processes (exponentially distributed power). 


2 

Var[v(tj^)J = (E V(t^))^ = V (t^^) 


(2.65 


Thus the tracking error variance is given by combining (2.63, (2,64, (2,65 
and (2,52 to yield: 


Va. 


E V (t.)(hj^-C g^) 


(2.66 
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For the model assumed in Fig. 2.7, 


V(t) - 1 


t< 


= 1+a 


t> 


^ T 


-T/2 -Mt 

T/2 = mt 
-T/2 < t < T/2 


(2.67 


For the square- non- overlapping gates assumed, let (At)'£-g, (At)Ug be the position 
of the lower and upper edge of the early gate, and let U^refer to the late 
gate. 


Then 


Z V(t.)h. 

1 L 



+ a/2 (-tg 




(2.68 


and 


S V(t.)g. = N N_ (1+a) 

I X p L 


(2.69 


where N = U -'t+lis the width of the early gate, is the width of the late 
E E E - ij 

gate, N is the number of pulses integrated, and V« is the mean power at the 
P S 

lower end of the early gate. 


The derivative of V(t) is given by: 

V(t) =0 t < -T/2, t > T/2 
|t|<T/2 


(2.70 
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Thus 


k = E V(t^)(h^ - C g^) 
i 


N E a/T 


p T 


p zmt 


Finally, (for non- over lapping gates) 


2 Ij 2 

Var e(0) = E V (t.) + 2 V (t^ 

r ^ 

E L 


But, V (t.) = 1 + a, ^ < i < U, , since the late gate is on the plateau. 

i 1* ■"* Ij 

Similarly, V(t^) = 1 + a/2 + (ia/2M)^ 5 ^ since the early gate is on 

the leading edge. 


Thus, 


var e(0) = N E (1 + | + ^ i) ^ Nj^N^d + a)^ 


= «p«e(\ + -k (k) 


2 (Ng-l)(2Ng-l) 


+ N^N (1 + a)^ 
L p 
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Substituting (2,62, (2,68 and (2.69 for C yields 


Var 6(0) = N^N ■ 2M 

^ Ju 


a- /as^ (N -1)(2N -dI 


+ - 4 - 


(2.74 


Substituting (2.74 and (2.71 ‘ into (2.52 yields the equation for the tracking 
error given as equation (2,53 in Section 2,3, 1.1, 
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2,3.2 


ContiTiiious Case 


In the previous section on the split-gate tracker, it was assumed 
that the data was sampled at time intervals far enough apart in range so that 
the sampled values were independent. In this section, the split-gate tracker 
equations are re-derived assuming that the sampled values are no longer inde- 
pendent. The continuous case (analog) is then approximated by letting the time 
samples occur at very closely spaced intervals. 


The derivation follows very closely the one given in the previous 
sections, thus only the modifications of that derivation will be given. 


In Section 2,3,1 it was shown that the variance of the tracking error 
for a single pulse is given by eq, (2.66. 


Var t “ Var e(0) 


(2.75 


where 


k = 9E e(t) 
9t 


t = 0 


(2.76 


is the mean error slope, and 


E(t^) = Z (2.77 

is the tracker error signal; 

Further, 

V(t^) is the value of the square law detected envelope at time tj^ 

h^ is the early gate response of the tracker 

g^ is the late gate response 
and C is a normalizing constant. 
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It was shown that k is given by: 



(2.78 


where T is the width of the leading edge of the return. This result is-not 
affected by correlation between the samples. 


The variance of e(0), however, must be recomputed. 

Now, 

Var e(0) = Z (h . ~Cg ) (h - Cg ) Coyar [V(t.)V(t )] (2.79 

ij ^ 

It can be shown, that for jointly Gatissian random variables, the covariance of 

the square law detected output is given by 

2 

Covar V(tpv(t^) = ^2.80 

where = E V(t^), and is the correlation between the Gaussian variables. 

If the early and late gates are separated enough to be independent, 
then the variance of e(0) becomes: 


Var e(0) 


I ViV p + C^(l+a)^ 
e.g. ^ 



(2.81 


where for simplicity , the gates are assumed to be square. The quantity, a, is the 
signal- to-noise ratio. The first and second sums above are over the early and 
late gates respectively. 

On the early gate, 


Let 



V. 

3. 

k 


= T 


( 6 + U.) 

1 


then 


(2.82 


(2.83 
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also 


C(l+a) = ^ 
k N, 


(2.84 


Thus, the tracking error becomes. 


Var t = 


Var e(0) 


= z (6 + u.)(6 + U.) 6^ E p? . 


e.g. 


j' lij 

N, 


N. 


6 T 


^ P. . + oS p . . 

e-8* xj ij 


N, 


N. 


+ 26T^ E n P, . + ^ U^U. P. 


i 11 

e.g»‘ Vr 

E 


e.g, 


ij 




(2.85 


If one assumes, stationarity , 


" Vj * 


( 2.86 


then it can be shown by symmetry arguments that ^ 


lUiPi^- =0 


(2.87 


Finally letting Z. = t./A where A is the width of the early gate, one has: 
lit t 


Var t 


n: 


E 

e.g. 


ij ^ 


+ 1 ^ P . . 

2 l.g. 


N, 




N, 


E Z Z p 

e.g.i J 


(2.88 
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This expression can be shown to be equivalent to that in Section 3.1 for 
independent samples. 

The summations in equation 2.88 have been evaluated for a 

Gaussian shape correlation function. These results are shown in Figure 2.9 

and 2.10. The results are given versus N the number of samples in the gate, 

and p the correlation between a sample in the middle of the gate, and one at 

1/2 

the edge of the gate. The results for N = 10, approximate fairly closely the 
continuous case. 

In the next section, the tracking error for the sampled split-gate, 
continuous split-gate, and MLE trackers are compared. 
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2.4 


Comparison of Tracker Performances 


In this section the accuracy expressions derived for the continuous 
split gate, the discrete split gate, and the MLE trackers are evaluated and 
compared . 

The continuous split-gate tracker accuracy expression of equation 
(2,88 requires a correlation function, P^j » between adjacent samples in the 
IF filter output. (It is assumed for simplicity that no pulse compression 
is performed) . The correlation function for the output process of the IF 
filter is given by the convolution of the correlation function of the input 
process with the autocorrelation function of the filter impulse response , xhe 
correlation function of the input process is just the autocorrelation function 
of the transmit pulse, since the return signal is the transmit pulse convolved 
with a spatially white distributed target. Thus if the IF filter is matched to 
a rectangular pulse of width 1/B, the output correlation function is the con- 
volution of four rectangular pulses. This function is approximately Gaussian 
in shape with variance given by 


= 2 


\l2 12 By 


( 2.89 


Thus, the IF filter output correlation, p(t), is approximately given by; 


p(x) = exp 


-3 (Bt) 


( 2.90 


[1+(BT)^] 
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Using this expression for P (T ) and BT=1 equation (2.88)is evaluated. The 
resulting performance is plotted in Figure 2.11 for several tracker con- 
figurations . 

The discrete split gate tracker variance is computed by substituting 

N = 1 into equation (2.88) and setting P . . = 0 for i ^ j. The performance of 
E 1 J 

this tracker is also presented in Figure 2.11. The discrete time model can be 
seen to be overly pessimistic relative to the continuous time model. 


A MLE tracker accuracy is obtained from equation (2.46) by 
substituting 


3 = 


2a, 


(2.91) 


and converting from time to distance units via 

j cAt 
d = 


(2.92) 


and 


a 

2 


Performing these substitutions yields for the MLE 

, 1/2 

1 — r“- r._ ' . 

a. 




(2.93) 


(2.94) 


This is also plotted in Figure 2.11 and compared with the split gate trackers 
(right ordinate) for a. = 2.5 ra and d - .5 m (T ~ 3 nsec). It should be noted 
that [k^ ^ 1 also depends on through (equation 2.48) 


o o 


U = O' 3 T 
max max 

= O' 3 At N 


(2.95) 


cells 


where = .3227 and number of resolution cells from the half 

power point of the return to the end of the observation interval. Substituting 
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-10 -5 0 5 10 15 20 25 


S/N in dB 

Figuie 2.11 . Comparison oi Tracker Accuracies 
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in (meters) for MLS and Split T^a-'keT-s (cL * 2.5 taeter, d = .5 meter) 




equation 2.91 and At = 3 nsec yields 

N 


U = .1452 -- 
max 


cells 


(2.96) 


- 1 ' 

In Figure 2.11, [k 1 is computed for N - 31 (corresponding to half 

t tz C S X X s 

o o 

the leading edge, 15 cells, plus 16 cells in the equivalent late gate) or 

U 2.0. A somewhat smaller tracking error can be achieved by the MLE 

max ~ 

by utilizing more cells on the plateau as illustrated in Figure 2.11 for the 
case =16. To compare the split-gate trackers at a particular waveheight 

with the MLE, use is made of the relationship between and T via 


a 



c 

2 


T 

3.1 


(2.97) 


which is obtained by a least squares linear fit to the leading edge of the error 
function. Thus converting the distance units, the error for the split-gate 
tracker is 


^ = 3.1 a^(NSD) (2.98) 

where NSD is the normalized standard deviation for the split-gate trackers 
(left ordinate) of Figure 2.11. 
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3.0 


BINARY PHASE CODE WITH DIGITAL PROCESSING 


3.1 SUMMARY 

The practicality of using phase-code and digital processing in a 
satellite altimeter and wave height analyzer is investigated in this rep,qrt. 

It is shown that a multi-bit digital processor (at least 2 bits plus 
sign) is required due to the distributive nature of the sea surface. 
Single-bit processing, which has been found effective in many different 

applications, is found to be unsuitable. Results show that for a compression 

* 

ratio of 1000:1 and an uncompressed return- to-noise ratio of -20 dB, the com- 
pressed output signal-to-noise ratio is only 4,1 dB which is insufficient to 
maintain the required tracking accuracy, 

3.2 PROBLEM DESCRIPTION 


A phase-code altimeter generates a phase-code sequence using feed- 
back shift- registers properly tapped as shown in Fig. 3.1. This sequence 
modulates the phase of a carrier. The resolution of the phase-code waveform 
depends on the duration of a bit, T, in the code and the output signal-to- 
noise ratio depends on the length of the code, K, as well as the reflected 
power-to-noise ratio. K is equal to the ratio of the duration of the entire 
waveform, T, to the resolution T i.e., 

T - Kt (3.1 

and it is often called the compression ratio. 

As defined here the return power is that from a single range cell* The 
equivalent peak-power- to-noise ratio is 10 dB for the example shown. 
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Fig. 3.1 Phase-code generator and transmitter. 


The receiver of a phase-code altimeter is shown in Fig. 3.2. As 


seen, the return 



Fig. 3*2 Digital phase- code receiving. 


is separated into two channels usually designated by I and Q for in-phase 
and quadrature. After mixing with the carriers, the video signals are 
quantized and delivered to a pair of digital-matched- filters , DMF* A DMF is 
a device which compares the input with the transmitted phase- code sequence 
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and totals the number of matches and the number of mismatches. The difference 
of the two numbers is used as an indication of the strength of the signal re- 
turn, Of course, in the absence of any return, this difference will be zero 
or nearly zero on the average, 

■-The^^NASA system requirements indicate that T is on the ‘ h- 
order of 3 nsec. Since the satellite is far from the target (556 km), the 
return is likely to be weak. The prsa-processed input signal-to-noise ratio 
for one range cell is expected to be on the order of -20 dB or less. Hence, 

K is necessarily long, say on the order of 1000. Together, the two requirements 
imply that the A/D must be fast enough to sample the input at 300 MElz rate and the 
DMF must be at least one thousand bits long. Obviously, if high accuracy 
quantization is also required, the hardware implementation of this receiver 
will be very difficult involving relatively high risk. 

In the following discussion, the performance degradation as compared 
to a linear system is analyzed as a function of the accuracy of quantization. 
For simplicity, the analysis has been limited to a smooth sea surface and a 
pulse-limited situation. The conclusions obtained and the method of analysis 
can be extended to other situations with a little modification. 

3 , 3 Basic Assumptions 

As indicated previously, it is assumed that the sea surface is flat. 

In such a case, the impulse response is a perfect step as shown in Fig. 3,3a. 
The return from a pulse-code waveform is a ramp as shown in Fig, 3^3b, com- 
posed of the returns from the different range cells. For a smooth sea, the 
power of the return, C, from the different range cells will be the same on the 
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Fig, 3.3a Impulse response of a smooth sea. 



Fig, 3,3b Response of an extended pulse waveform. 


average. The receiver, in effect, cross correlates the entire return with a 
replica of the transmitted pulse-code sequence once every t seconds. The in- 
stant such a sequence in phase with the replica appears at the return, the 
cross correlator output will peak indicating its presence. Since the return 
is an overlap of many pulse-code sequences with different phase shifts, the 
cross correlator output will also have smaller peaks due to partial matching 
between the replica and the return# The power of the subpeaks Is proportional 
to the sidelobes of the auto-correlation function of the phase-code sequence. A 
t 5 rpical auto-correlation function of such a sequence is shown in Fig, 3,4. 
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These subpeaks are undesirable and they may be reduced by proper signal 
design amplitude tapering of the replica. In the present analysis* it 
is convenient to assume that the sidelobes are identically zero. To a 
first approximation, this assumption does not affect the conclusions to 
any large extent* 

It is assumed that the return from a range cell has the form c s(t) 
where c is a complex Gaussian process and s(t) is the transmitted signal, 
c is independent from range cell to range cell and the amplitude of c is 
Rayleigh whose mean power is C, Due to this asstimption, the compressor output 
signal" to~noise ratio of either branch will be the same as that of the entire 
processor as shown in Fig. 3* 2, 

It will also be convenient at this point to define the s)rmbols to be 
used with the A/D converters. The converters considered here are always odd, 
uniform, and without any dead zone. The input-output relationship of such 
a quantizer is shown in Fig, 3.5, 

Tapering will be more difficult in digital systems especially in single-bit 
systems. It is expected that larger signal power loss will result in sup- 
pressing sidelobes in digital systems. 
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Fig. 3.5 A quantizer. 


Shown is a two-bit (plus sign) quantizer with M = 4 output levels on the 
positive input side. The spacing between two output levels is A = 1.0 
and the spacing between any two adjacent input thresholds is ”a". 
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3.4 


Digital Compressor Output Power 


As a measure of the system performance, the output power of the digital 
processor (one branch) will be computed and compared with the output power of a 
linear processor using phase-codes. Intuitively, it is expected that a digital 
or linear compressor restores the return from that shown in Fig. 3.3b to that shown 
in Fig, 3,3a, The sharpness of the edge and the height of the jump of the recon- 
structed signal is an indication of the effectiveness of the compressor. As will 
be shown later, there is a loss of the height of the jump as the quantization 
error increases. 

Let s^, i - 1 ... K, denote the phase-code sequence, s^ may be plus 
or minus one. Let v^ denote the video input sequence^w^ denote the video quantized 
sequence and y^ denote the compressed output sequence. Assume the time origin 
is centered on the instant prior to the arrival of the earliest echo. Then, the 
video input sequence v, is given by 


^ K(j) 

V, - n. + 2 

J j i=i 




(3.2 


where n. 

J 


is the Gaussian noise, 


c . 
1 


is the reflective coefficients and 


K(j) =0 j < 1 

= j K > j > 1 

= K j > K 


The compressor output is given by 

j+K 


K 


V . . = 2 s w. 

n=l ^ 


(3.3 • 
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since both n. and c. 

3 3 

w. also has mean zero 


3 


are Gaussian of mean zero, and the quantizer is odd, then 
, The compressor output power is equal to the expected value 
this power is 




K 

S E C(s 

n 

n,m 


"jW <*m 


j+m 


)] 


K 

Z E 
n-1 


j+n 


] + 


K 

E E 
m?^n 


j+n 


w. , ]s 
j+m n 


m 


(3.4 


The expression on the RHS is split into two sums because it is easier to handle 
them separately. The terms in the first sum are simply the average quantizer 
output powers. Each of the K terms in this sum may be computed by the basic 
relationship 


03 

e[w^J = J p.(w ) 
3 ^00 3 


dw 


(3.5 


where p.(w) is the probability density function of w.. In the present case, w. 
3 3 3 

Is discrete so that the integral in (3.5 can be replaced hy a finite sum as 


2 M M 

E [w ] = 2 l; p (L ) = 2 S l; p (L.) (3.6 

3 i=-M -^3 1 ^ ^ 


where each L. is a quantizer output level (see Fig. 3,5 and the probability p.(L.) 

^ 3 i 

is given by the probability that the quantizer input v^ falls between the i 


and (i-1) 


threshold. Since v^ as given in (3,2 is Gaussian of mean zero. 
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p ap = 0.5 - Q(a/Oj) 


= Q ((i-l)a/aj) - Q (ia/a.) 

1 < i < M 


Pj<V = Q ((M-l)a/aj) (3.7 


2 

where is the variance of (see (3*11 below). The function Q (x) is 
the tail of the normal probability intergral, and is given by: 


Q (X) = 




QQ 2 


(3.8 


The terms in the second sum in (3.4 are the cross correlations between 
pairs of quantized Gaussian random variables. They are related to the input 
correlations as 


R (w w ) 

m n 


M-1 

4 *?] E fA exp 
1^1 u=*o' 

I odd 



X I A exp 
u 



A = 0.5 

o 


R (v V ) 
m n 

I I 

2tt a_ a„ 't ! 
m n 


(3.9 


A = 1, u > 1 

u ' — 
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where H (x) is the Hermite pol 3 momial of ('f'-l)th order. This relationship 
'C/*' 1 

can be derived from the general result given in [ 1 ] using the same approach 
as in [2]. The details are omitted in this report. 


At this point, (3,9 can be substituted directly into (3.4 to find 
the output power of the digital compressor but considerable simplification 
is possible due to the fact that R(v^ ^n small, for m n. From (3,2 

m n 


K(m) K(n) 

R(v V ) E [n n ] + eQ Sc s. e - - s, ] 

^mn ‘'rfc m-i+1 i n-k+1 k 


K(m) K(n) 

= N 6^ + S S E ^n-k+1^ 


®k 


(3.10 


When n = m, R(v , v ) = a , the variance of v , and is given by: 
’ n m m’ m' o ^ 


K(m) K(m) 


a = N + 2 

m i 


^ '- 1.1 ^ i_i_i = N + K (m) C 

m-i+1 n-k+1 


(3.11 


2 

using the fact that E[c, c.] = 0, i ^ j and e[c. ] = C. On the other hand, when 

i J ^ 


m 5 ^ n, R(v v ) is 
’ ' m n 


R(v 


m 


v ) ^ 
n 


K' (m,n) 

C 2: s 

k==l 


m-n+k k 


where K'(m,n) = min (K(m), K(n)) 


(3.12 
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For a well designed code with small sidelobes, the result of the summation in 

(3.12 will be small but in any case it is always less than K(n). Hence, 

R(v V )/ct c is small if K(n) » 1 or if N » K(n) G. In practice* one of the 
in n 111 n 

two conditions is always satisfied. 

Since (3.9 is a power series in R(v v )/a a . the first order term 

m n m n 

in (3,9 will be a good approximation for m ^ n. Substituting (3,10 and (3.9 :> 
into (3.4 it is found that 




K 

E 

n=l 


+ E 


E s s 
/ n m 
mfn 




where 


K(j-Hn) K(j+n) 

E E 

i k 


^j-hn-i+1 '^j+n-k+l ^i ®k ' 


(3.13 


M-1 / 2 2\ 

f(CT)=2 E A exp ^ » ) / /5 tt a 
n=0 ^ \ 2o / 


Adding and subtracting the term E 
(4,12) can be put into the form 


K ^ « K ( j+n) 

"j+n-i+1 "j+n-k+1 ^1 \ 

n— 1 1 , K 


= S E [w ‘ 
n=l 


K 


■] - S 

n-1 




K(j+n) 

E C 


/ K 
+ E E 
\ n-1 


s 

n 




K(j+n) 

E 

i=l 


^j+n-i+1 



(3.14 


The first two sums in (3,14 have the following interpretations. The first term 
is the total quantizer output power in an interval T = Kt . The second term is 
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the output target return power over the same interval. The difference of the 
two terms is, therefore, the total output noise power. 

The third sum, denote it by B, can be simplified considerably. Letting 
T = n-i, it becomes 


B = e( 2 2 


T n=l 


The limits on T depends on j as follows. 


(3.15 


T = -(j-1) ,K+j-l 


j < 0 


j , . • • , K+ j j ^ 0 

K 

The summations S f(a,. ) s s are the sidelobes of the quantized signal 
- 1 +n n n- t -i » 

n=l 

when T 0, As a result of the zero sidelobe assumption, B is reduced to 


B 


K 

S c . 
1 J+1 
,n=l 



2 


^ C 



j+n 



(3.16 


for j ^ 0 and is zero for j < 0, Thus, the third sum in (3,14 is a measure of 

the compressor output due to the return which is in-phase with the local signal-replica. 

In the next section, the response of a hard!- limiting (one bit I and 
one bit Q, including sign) system will be analyzed in detail to be followed by 
the analysis of multi-bit system performances. 
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3.5 


Hardliiniting Digital Processor 


Since the statistics of the input to the receiver is non-stationary, 
it will be convenient to compute the compressor output for different time in- 
tervals separately. Again, as above, j-1 corresponds to the instant prior to the 
earliest return reaching the compressor. Fig. 3,6 shows the breakdown of the 
time segments. 



Fig, 3.6 Breakdown of time segments. 

Case (!) j+K < 0 

Prior to the instant j= -K, the input in the compressor is en- 
tirely Gaussian noise generated by the front-end of the receiver. Since this 
noise is uncorrelated with the transmitted signal, the second and third term in 
(3,6 sre zero. The compressor output is given by the first term only, which 
may be computed using (3,6, (3,7 and (2.8. For a hardlimiting system, this 
is quite simple since all positive input is mapped into (+0,5) and all negative 
input is mapped into (-0,5), The output power is, on the average 
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(3.17 


= K(0.5)^ = K/4, 


j < -K 


For a hardlimiter, the first term in (3.14 is always K/4. 

Case (ii) 0 < j+K <K-1 

( 

During this time interval, the target return is growing as j increases 
and is competing with the Gaussian noise for quantizer power. There is as yet 
no coherent output power at the compressor and the third term in (3.14 is still 
zero. The terms in the second sum are, however, growing linearly as a function 
of j. The average compressor output is 


9 j+K 9 

E ] = K/4 - Z i C f^ (a,) 


j+K 


(3.18 


i=l 


where some simplification has been made based on (3.14 and the definition of 
the upper limit K(j+n). is given by (3.11 as 


= i C + N 

1 


(3.19 


f(a^) is 


1 1 
" /2tt a. "" /27T (~i”c + N) 


(3.20 
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Consequently, 


® = K/4 


2 

1=1 


21 C 

n(l C+ N) 


(3.21 


Case (111) K < j+K < 2K-1 

2 

From the time j-1 on,E[j^^,+K] will always be a return component which 
is matched to the transmitted-signal-replica. However, the share of the quantizer 
power devoted to this particular component is gradually changing from a peak at 
j=l to a constant level at j=K, The total input power is 


N + j C 


j < K 


^ N + K C 


j > K 


(3.22 


so that 


£(a^) 


1 

/2TT(j C+N) 


j < K-1 


The output compressor power is 


(3.23 


E [y ] = K/4 - i 


s — Li.g_ ]+ 2 k c (1-1) 

i=j n(jC + N)/ tt(KC + N) 


+ \ 


2(j-l) VF ^ I IJZ 


/2tt(K C + N) y5n(i C + N) 


(3.24 
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Case (iv) j + K > 2K-1 


Beyond j “ K, the return from the sea surface has reached a maximum 


and the average power will remain constant. ^\y ^.v ^ is simply 

j+K 


E Cy ] = K/4 - % 


2K K C 
tt(K C + N) 


+ % 


2k/C 
v^^'(K C + N) 


= K/4 


(3.25 


Let K = 1000. The digital compressor output as a function of time is shown in 
Fig. 3.7. 


B 


250 



4 

-1000 




250 


■i Time 

1000 


Fig* 3.7 Average compressor output. 

The power levels reached in Fig. 3,7 depend on the ratio C/N, For different 
values of (C/N) these levels are tabulated in Table I, 
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TABLE I, Average Compressor Output Levels 


(C/N) 


A 


Coherent 
Power (B~A) 


10 

98 

616 

518 

30 

108 

557 

449 

50 

115 

522 

407 

100 

129 

470 

341 

200 

147 

415 

268 

300 

160 

384 

224 

500 

178 

348 

170 

1000 

201 

310 

109 


For a given K, a linear compressor output will be 


E [y ^ KN 


j < 1 


= KN + K C j ^ 1 


(3.26 


To compare various system performances, it will be convenient to define a 
'jump' ratio as 


Y = 


E - E [y/] 

E [y^] 


(3.27 
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For the linear case, this ratio is simply 


linear 



(3.28 


For the hardlimiting case, tabulated in Table II for various values of 

N/C with K ■= 1000 


TABLE II. Y, . And Yut K. = 1000 

j-ij.n6d.ir riJ-i 


(C/N)'^ 

Y 

Linear 

% 

Loss (dB) 

10 

100.0 

5.29 

12.77 

30 

33.3 

4.18 

8.95 

50 

20.0 

3.54 

7.52 

100 

10.0 

2.65 

5.77 

200 

5.0 

1.81 

4.41 

300 

3.3 

1.38 

3.77 

500 

2.0 

.96 

3.18 


The degradation defined as Yjj^ plotted in Fig. 3,8 

as a function of N/C. It should be noted that this loss in performance cannot 

be compensated for by raising the transmitter power by a factor equal to the ratio 

\inear^"*^HL* seen, even though a system with N/C = 100 and 

K = 1000 has a loss of 5.77 dB, an N/C ten times smaller can only bring the hard- 

limiting system performance to about half that of the linear system. In fact, 

for large K it may be shown that Yur is approximately 

HL 
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2 


(3,29 


'HL 


— h h 

TT [(K+B)^ - (1 + P)^] 

K(1 . ^ + 3(^n(K + 3) - ^n(l + 3)) 


K » 1 
g = N/C 


Therefore, even if p = 0 


Y„ = - / (1 - |) = 7.07 3 -* 0 

'HL TT TT 


(3,30 


Thus, the degradation as compared to a linear system is inversely pro- 
portional to 3 as g 0, On the other hand, (3.29 also shows that as 

3 '”» ^HL 


= 1 

m TT 3 


3 ® 


so that the degradation approaches — (2 dB) . 

TT 


' (3.31 


In view of the present NASA requirements, y must be close to 
10 dB to achieve the necessary tracking accuracy, A hardlimiting system is, 
therefore, inadequate. In the next section, multi-bit systems will be con- 
sidered. 
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3.6 


Multibit Digital Processor 


Since in the present application the compressor output at the 'jump' is 
the most important, the ratio Y 3,s defined in Section 3.5 is the only quantity 
being considered in detail here. Generally, the compressor output will re- 
semble that shown in Fig, 3,7 for 1 or 2-bit systems and will approach the ideal 
step function as more bits are added. 


The equations needed are (3,14 and (3,16, They are evaluated 

numerically on a UNIVAC 1108 computer. The results are shown in Figs. 3,9 

* 

and 3,10 for one to five-bit processors. The K and N/C assumed in these com- 
putations are 1000 and 100 respectively so that a linear system will have a 
''^Linear these figures, y is plotted as a function of the normalized 

threshold spacing a//C, As seen, there is an optimum spacing for each of the 
five processors. Generally, the optimum is quite broad. The best that an M- 
bit system can obtain is shown in Table III,, 


TABLE III, Digital Processor Performance, y (K ~ 1000, N/C = 100) 


HL 

2.65 

1-bit 

5.40 

2-bit 

7,85 

3-bit 

9.11 

4-bit 

9.70 

5-bit 

9.92 

Linear 

10,0 


Sign bit is extra. 
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3.7 


Beam Limited Case 


It was shown in the previous sections that there 
is a significant loss (6 dB) associated with hard limiting processing. 

The model used to arrive at this result assumed that the altimeter - uncompressed 
return is pulse limited. In practice, the occurrence of this situation depends 
on the choices of pulse duration and antenna beam width. Very often, the system 
parameters are so chosen that the return is only partially pulse-limited. As a 
result, the system loss will depend on the pulse limiting factor P defined on the 
ratio of the number of range cells in a single pulse versus the number of range 
cells being illuminated by the antenna at the sea surface. For compression ratio 
K = 1000 and noise- to-slngle cell return of 100, the loss as a function of p is 
plotted in Fig 3.11. As seen, the loss is above 5 dB when 3 is larger than .3. 
Thus, this new result does not affect previous conclusions significantly. 



Pulse Limiting Factor 
K = 1000, N/C = 100 

Fig, 3,11 Loss as a function of pulse limiting factor. 
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Based on these results a 2-bit (plus sign) digital processor seems adequate for 
the present application and the optimum threshold spacing is a//C = 15. The 
compressor output power for the entire time interval of interest is shown in 
Fig. 3.12. 


525 


4633 



1 


4633 


Time 


Fig. 3,12 Compressor output of a 2-bit (plus sign) digital processor 

(K = 1000). 


3.8 Conclusion 

Analytical results have been obtained in this report to show that a 
hardlimiting digital processor operating on a phase-code signal is inadequate in 
the NASA altimeter and wave height measurement applications. It is shown 
that because of the distributive nature of the sea surface, the 

hardlimiting processor output signal- to-noise is bounded at a low value regardless 
of the compression ratio K or the return echo-power-to-noise ratio as long as K is 
much larger than unity. In the example shown, the output signal- to-noise ratio is 
only 4.1 dB for K = 1000 and C/N = -20 dB. This is insufficient to maintain the 
required tracking accuracy. 

It is also shown in this report that the digital processor performance 
improves as the quantization error is reduced, A 2-bit-plus- sign processor 
seems to be adequate when the threshold spacings are properly set. However, 
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at the speed of the intended operation and in a satellite environment, there 
is much risk involved in implementing such a complicated digital processor. 

The cost, as compared to a linear FM system, is also likely to be much greater. 
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4.0 


COMPARISON OF SEA SURTACE CROSS SECTION MEASUREMENTS AT NORMAL 


INCIDENCE WITH BARRICK'S MODEL 

In this appendix, Barrick's model for the normalized sea cross-section 
is compared with several sets of experimental results. It is found that his 
model predicts values for the cross-section per unit area which are somewhat 
higher than the measured values and suggests that reducing Barrick's model by 
about 6 dB would give a more conservative model. 


Barrick [1], [2] gives the cross section per unit area of the sea 
as : / 


where 


and 


0 4 / 

^ 2 ” 

0 = Sec Y / exp 

-tan Y 

2 1 

^2 

s \ 

L " J 


(4.1 


/, 


Y is the angle of incidence (i.e. angle from vertical) 


S is the nns sea slope* 

Barrick relates the rms sea slope to wind velocity and wave height by: 

= 5.5xl0“\ and = 2.55 
h 


(4.2 


where V is in m/sec. Barrick's model is based on specular scattering, and 
thus has no frequency dependence. At normal incidence e.g.(4.1 gives. 


a‘^(dB) = 22.6 - 10 log V(m/sec) 

= 19.4 - 10 log V(kt) 


(4.3 


or 


a°(dB) = 13,6 - 5 log o (m) 


(4.4 
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In Figs. 4,1 - 4.6, equation (4.3 is compared with experimental measure- 
ments as reported by several references. In Figs. 4.2 - 4.4, the median cross- 
section was reported. For those figures, eq. (4.3 is reduced by 1.59 dB, This 
corresponds to the mean-to-median ratio for Rayleigh clutter. 

Three comments seem appropriate to Figs. 4,1 - 4.6. 

a) There is a large scatter to the experimental data which is 
probably not entirely due to experimental error, (e.g. re- 
duction of the measured values due to finite beam width) . 

b) There appears to be a trend towards reduced cross-section 
at higher wind speeds. Further, Barrick's model reflects 
this trend moderately well. 

c) Barrick's model seems to be somewhat optimistic, particularly 
when compared to the NRL data on Figs, 4,2 - 4.4. 

These comments suggest that Barrick's model 4.3 can be taken as a 
relative optimistic prediction of the mean cross-section at normal incidence. 

a^(dB) = 16,6-10 log V (m/sec) 

= 7.6-5 log (m) 

might represent a more conservative estimate. 

The scatter of the data strongly suggests that any proposed altimeter 
design should not be severely degraded by even a 10 to 15 dB reduction of 
relative to Barrick's model. 
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O Campbell 3,4 cm [6] 
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Georgia Inst. [12] 


O 35 GHz 

® 24 GHz 
^ 8.8-10 GHz 





4.1 


Effect oti Tracker Performance 


To determine the effect of sea state and the model used on tracking 
accuracy, the performance of several split-gate trackers was computed. The 
results are given versus sea state for Barrick's model, and for Barrick's 
model minus 6 dB, To make the comparison, a signal-to-noise ratio at some 
value of a° (cross-section) had to be assumed, A design point of S/N = 10 dB 
at = 6 dB was used. The results are given in Figure 4,7, 

As was expected, the tracking accuracy is affected by the reflectivity 
model used, however, the effect is not severe. There is about. 30% degradation 
between results using Barrick's model, and results using the more conservative 
model. As a point of reference, for a significant wave height of 10 m, Barrick's 
model gives == 11.6 dB, while the conservative estimate gives = 5,6 dB, 
Recall that the design value is = 6 dB at a 10 meter ;^ave height; thus the 

conservative model is being used in the altimeter design. 

It is interesting to note that if the tracking accuracy is normalized 
to significant wave height, the results are almost independent of wave height 

above 3 m. Thus the signal loss at the higher sea states is not significantly 

affecting the tracking accuracy. The most important effect is simply that the 
leading edge of the return becomes less distinct at high sea states. 
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FIGURE 4.7 

Comparison of Tracker Performance versus Wave Height for 
Two Sea Reflectivity Models* 


These numbers are normalized to significant wave height rather than 
mp width, 1.29 (CT/2)). Further, they do not include 

the effect of correlation between samples, 
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5.0 


AIRCRAFT SYSTEM PARAMETERS 


Since the system design was based on orbital parameters and whereas the 
breadboard unit will actually be tested using an aircraft, it is necessary to 
modify the selected set of system parameters given in Part I. This section describes 
the required changes in the system parameters for an aircraft mounted altimeter. 

It is shown that the reduction in height and velocity from a .satellite to an air- 
craft requires an increase in antenna beamwidth from 3^ to 15° and a reduction in 
peak power from 2 kw to 1.0 watts. A reduction in the maximum spatial frequency 
(associated with the Gulf Stream) to ,007 Hz enables the tracker bandwidth to be 

reduced to about .06 Hz. In order to check the validity of the predicted altimeter 

performance (and for investigative purposes) it is suggested that some of the 

system parameters be made variable over a range of values. Specifically, a recom- 

mended range of values are given for tracking gate widths, compression ratio and 
pulse length, PRF and bandwidth. 


- 172 - 



It would be desirable to study, with the aircraft-mounted altimeter, 


the altitude measurement accuracy and how it is affected by changes in: 

1. Tracking gate widths; 

Early gate: 1-8 resolution cells 
Late gate: 1-16 resolution cells 

2. Compressed pulse length: 3.0 nsec to 10 nsec 

3. PRF: vary from half to twice the inverse decorrelation time 

4. Number of uncorrelated pulses averaged by the loop filter 

These test objectives require modifications to the orbital values of 
system parameters as shown in Table 5.0-1 and described in the following sections. 


5.1 Beam Width and Tracker Gate Widths 


A minimum of 24 resolution cells must be contained within the return 
beam-limited pulse length, t^, in order to test the suggested ranges of gate 
widths. At the maximum compressed pulse width of 10 nsec, this requires a 240 
nsec surface return pulse lengthy This pulse length is given by 


t 


B 




where h is altitude, c is speed of light, and 6^ is 3 dB beamwidth of the altimeter 
dish. At an aircraft altitude of 3048 m, a beamwidth of 17.5° is necessary for 
a 240 nanosecond t_. This is in contrast to the satellite case, in which a 3° 

D 

beamwidth at an altitude of 556 km provides a t^ of 1,2 |isec. Thus the 
antenna beam must be opened up for the aircraft altimeter, in order to approximate 
the duration of the satellite sea surface return. 
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Table 5.0-1 System Parameters (Aircraft Mounted Altimeter) 


I. 

Radar Parameters 



Antenna Beamwidth 

15^ 


Pointing Accuracy 

2.5° 


Antenna Gain Peak 

Average 

20.9 dB 
20.25 dB 


Peak Power 

200 raw (min) 

(use 1 watt TWT) 


System losses (other than processing 
losses in pulse compressor) 

5 dB 


Noise Figure 

5.5 dB 


Frequency 

13.9 GHz 


Pulse Compression Processing Loss 

.55 dB 


*Main Lobe Broadening Due to Tapering 

23% 


* 

Compressed Pulse Length (after tapering) 

3 to 10 nsec 


Tracker Bandwidth 

> .06 Hz 


PRF 

100 to 1000 Hz 


Number of Samples to be Integrated 
by Tracker 

700 to 3400 


Single Pulse SNR 

10 dB 


Compression Ratio 

irk 

1000 to 100 


Bandwidth 

kk 

360 MHz to 120 MHz 


Uncompressed Pulse Length 

** 

2.0 jjtsec to ,9 p.sec 

II. 

Ocean Cross-Section 

= +6 dB 

III. 

Maximum Spatial Frequency 

SF = .007 Hz 

IV. 

Aircraft Parameters 



Height 

3.0 km 


Velocity 

71 m/sec (140 knots) 


Modified Taylor weighting, -25 dB peak sidelobe, to be switched in or out. 

Variable compression ratio, bandwidth, and uncompressed pulse length to be 
obtained by generating maximum linear FM waveform and gating out desired 
segment. 
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It should be noted that, at the lowest sea states, the surface shaping 
function will limit the effective beamwidth to less than 17,5°, so that the full 
24 resolution cells will not be available in the aircraft test for all sea states. 
Consequently, a 15° beam is suggested for use in the aircraft version. The surface 
shaping function limitations are summarized in Table 5.1-1- 


TABLE 5.1-1 

SURFACE SHAPING FUNCTION T, IMITATIONS 


Wind Speed 
(knots) 

Effective 3dB 
Beaim-7idth 

Surface-Shaping-Limited Pulse 
Duration at 10,000' (nsec) 

7 

13,5° 

140 

10 

16.1° 

200 

14 

19° 

280 

19 

22° 

380 

22 

24° 

450 


5,2 Compressed Pulse Length 


For a surface wave device pulse compressor, the compressed pulse 
length could be varied by interchanging surface wave devices of differing band- 
widths, A simpler and more economical technique, for a linear FM waveform, would 
be simply to gate out a portion of the transmit pulse, thereby decreasing the 
dispersion and bandwidth proportionally. This technique has the further advantage 
of keeping the single pulse signal-to-noise ratio at the pulse compressor output 
constant as bandwidth is changed. The required compression ratio (CR reg’d) for 
a given signal-to-noise ratio is 


V 7 

CR . , = = KB 

req*d 2 
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where t is compressed pulse length, B is waveform bandwidth. Thus, as B is 
c 

2 

varied, a constant signal-to-noise ratio is maintained if CR is varied as B , 
This variation of CR is realized if the dispersion (uncompressed pulse length) 
is varied proportional to B, since CR = BT, 


5.3 PRF 

As a result of the decrease in altimeter velocity, the decorrelation 

time of the surface return has increased significantly. For a 3.0 nsec compressed 

* 

pulse, the decorrelation time for a 140 knot aircraft is 4,8 msec , compared to 
,6 msec for the same pulse length in a satellite. Thus the decorrelation- limited 
PRF is about 200 Hz at this pulse length whereas at a 10 nsec pulse length, the 
decorrelation- limited PRF is 360 Hz. The PRF should therefore be variable from 
100 to 1000 Hz, to cover a range of half to twice the inverse decorrelation time 
at all pulse lengths. 

5.4 Pulse Averaged by Loop Filter 

The number of pulses averaged by the tracking loop filter is roughly 
the ratio of the PRF to the loop bandwidth. Thus, provision should be made to 
vary this ratio, as well as the PRF, 

5.5 Peak Transmitter Power 

Even with a 15^ bearawidth, the power in the surface return would be 
40 dB higher than that for the satellite altimeter with a 3° beam. This is 
because of the large difference in altitudes, since the return signal strength 
from the pulse-limited footprint varies as one over the cube of the altitude. 
Consequently, a savings can be realized in the aircraft breadboard model by elimin- 
ating the final TWT in the transmitter and using a 1 watt TWT, 

* 

Wave spray may decrease this number to about 1,8 msec. 
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